Accepted to The Astrophysical Journal on June 9, 1998 



AN OPTICAL/NEAR-INFRARED STUDY OF RADIO-LOUD QUASAR 
ENVIRONMENTS II. IMAGING RESULTS 

Patrick B. Hall Q | 
Steward Observatory, The University of Arizona, Tucson, Arizona 85721 
Electronic Mail: hall@astro.utoronto.ca 

Richard F. Green 

National Optical Astronomy Observatories, Tucson, Arizona 85726-6732 
Electronic Mail: rgreen@noao.edu 

ABSTRACT 

We have previously reported a significant excess of K^19 galaxies in the fields 



of a sample of 31 z=l-2 quasars ( [Hall, Green fc Cohen 1998; ), Here we examine the 



properties of this excess galaxy population using optical and near-IR imaging. 

The excess occurs on two spatial scales. One component lies at #<40" from the 
quasars and is significant compared to the galaxy surface density at 6>40" in the same 
fields. The other component appears roughly uniform to #~100" and is significant 
compared to the galaxy surface density seen in random-field surveys in the literature. 

The r—K color distributions of the excess galaxy populations are indistinguishable, 
and are significantly redder than the color distribution of the field population. The 
excess galaxy population is consistent with being predominantly early-type galaxies at 
the quasar redshifts, while there is no evidence that it is associated with intervening 
Mg II absorption systems. The average excess within Q.bh^ Mpc (~65") of the quasars 
corresponds to Abell richness class ~0 compared to the galaxy surface density at 
>0.5/iy 5 1 Mpc from the quasars, and to Abell richness class ~1.5 compared to that 
from the literature. We estimate -O.65I055 magnitudes of evolution in to z=1.67 
by assuming the excess galaxies are at the quasar redshifts. 

We discuss the spectral energy distributions (SEDs) of galaxies in fields with data 
in several passbands. Most candidate quasar-associated galaxies are consistent with 
being 2-3 Gyr old early-types at the quasar redshifts of z~1.5. However, some objects 
have SEDs consistent with being 4-5 Gyr old at z~1.5, and a number of others are 
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consistent with ~2 Gyr old but dust-reddened galaxies at the quasar redshifts. These 
potentially different galaxy types suggest there may be considerable dispersion in the 
properties of early-type cluster galaxies at z~1.5. There is also a population of galaxies 
whose SEDs are best modelled by background galaxies at z^>2.5. 



Subject headings: Surveys — Quasars: General — Galaxies: General, Clusters of 
Galaxies 



1. Introduction 

With continuing advances in instrumentation now enabling detailed studies of very faint and 
distant galaxies, it is useful to seek efficient methods to find galaxies and clusters at z>l. One 
such possible method is to look for galaxies associated with quasars, specifically the radio-loud 
quasars (RLQs) which comprise <^10% of the quasar population. 

Radio-quiet quasars (RQQs) are rarely found in clusters at any redshift, but ~35% of 
intrinsically luminous (Mb<— 25) RLQs are located in clusters of Abell richness class 0-1 (and 
occasionally 2) at z=0.5-0.7 ( |Yee &: Green 1987 ). These quasar host clusters typically have 



anomalously low X-ray luminosities Lx (Hall et al. 1995; Hall, Ellingson & Green 1997) and 



velocity dispersions a v for their richnesses, and thus may be younger and less virialized than 
optically-selected clusters flEllingson, Green fc Yee 1991 ; hereafter EGY91). In addition, RLQ 



environments are known to evolve rapidly and differentially. At z<0.5 only low luminosity RLQs 
are seen in richness 1 clusters, whereas at ,2=0.5-0. 7 both high and low luminosity RLQs can 
be found in such environments. A similar effect is seen for FR II radio galaxies ( Hill fc Lilly| 



1991). EGY91 model this evolution by postulating that galaxy interactions (more frequent in 
young, low-<7 v clusters) are largely responsible for creating and fueling RLQs, which then fade as 
the interaction rate decreases in their evolving, virializing host clusters. Yee & Ellingson (1993) 
suggested that if the cluster formation rate dropped at z~0.7, most quasars we see in rich clusters 
at z<0.7 would be old ones fading on the host cluster's dynamical timescale, and that luminous 
RLQs might be found in rich clusters at z>0.7. The outstanding feature of quasar evolution is the 
sharp peak in space density and/or luminosity at z~2-3 ( |Shaver et al. 19*96 ), but little work has 
previously been done on quasar environments at z>0.7. 

Some RLQs show possible additional evidence for being located in rich environments, in the 
form of an excess number of "associated" C iv ( |Foltz et al. 198S| ) or Mg n ( [Alder oft, Bechtold & 



Elvis 1994 ) absorption systems within ±5000 km s 1 of the quasar redshift. There is a tendency 



for associated absorption to be preferentially found in steep spectrum sources. Anderson et al. 
(1987) quote a 2a preference for associated C IV in the unpublished Radio-Loud Survey of Foltz 
etal. This can be seen in Fig. 3 of Foltz etal. (1988): 16 of 22 RLQs with strong associated 
C IV absorption have steep radio spectra, whereas only ~10 would be expected. As discussed in 
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Aldcroft, Bechtold & Elvis (1994), this dependence is also suggested by the studies of Foltz etal. 
(1986) and Sargent, Boksenberg & Steidel (1988). Associated absorption systems may arise in 
gas expelled at high velocity from the quasars or in galaxies in clusters at or near the quasar 
redshifts. In the last few years detailed spectroscopy has shown that some associated C IV systems 
are almost certainly intrinsic to the quasars, as discussed in the introduction to Paper 1 ( [Hall 
Green fc Cohen 1998| ). However, it remains possible that a substantial fraction of such quasars 
reside in clusters which produce associated absorption. Spectroscopic and imaging approaches are 
complementary ways to investigate this question. 

In 1994 we embarked upon a project to study the environments of RLQs to z=2.0, to study 
any correlations between RLQ environment and quasar properties such as associated absorption, 
and to study any examples of high-redshift galaxies and/or clusters found in RLQ fields. Paper 1 
( Hall, Green Sz Cohen 1998| ) describes our sample and observations as well as our data reduction 
and cataloguing techniques. We observed <^3'x3' fields around 33 RLQs with z=l-2 and two 
control fields with typical 3a limits of r=25.5, i^ s =20.5 for z<lA, and K s =21 for z>1.4. An 
unusual aspect of our data reduction procedure is the creation of images with constant RMS 
noise from mosaiced images with exposure times that vary from pixel to pixel. A 3a detection 
thus represents a brighter galaxy at the edge of an image than in the center, but as long as 
the exposure time is everywhere >36% of the maximum, the average 5a magnitude limit will 
be greater than the 3a detection limit across the entire field. Object detection was done on 
RMS-weighted sums of the r, K s , and J data (where available). Photometry was done in all 
filters using the detection aperture from the summed image. We use FOCAS total magnitudes 
and colors determined from FOCAS isophotal magnitudes. The useful area surveyed consists of 
221 arcmin 2 of RLQ fields and 19 arcmin 2 of control fields. Our control field galaxy counts (and 
color distribution; see Appendix [A]) are consistent with the literature average, but the statistical 
uncertainties are unacceptably large for our purposes, so we also assembled a comprehensive 
set of random-field galaxy surveys from the literature for comparison. One complication arises 
from comparison with this literature sample: our 1994 and 1995 photometric solutions produce 
a -fT-band galaxy number-magnitude relation N(m) slightly higher than the literature at bright 
magnitudes. Given the possible systematic differences between the various iT-band surveys, we 
conduct some of our analyses using a conservative magnitude scale (defined so that our number 
counts exactly match the literature around K=17) as well as our original liberal magnitude scale. 
Both scales use UKIRT K magnitudes; the only differences between them are systematic offsets: 
on the conservative scale the 1994 magnitudes (mostly z>lA objects) are 0^06 fainter than the 
liberal scale, and the 1995 magnitudes (mostly z<1.4) are CP12 fainter, for an average difference 
of ~0? 1 08 given the numbers of galaxies in each subsample. As seen in Figure |l|, even under our 
conservative magnitude scaling the galaxy counts at K s ^19 in 31 RLQ fields with good data are 
higher than those of field surveys; i.e., there is an excess of galaxies in our combined RLQ fields. 

This paper analyzes the excess galaxy population found in Paper 1. In §|| we show that there 
is an excess of galaxies spatially concentrated around the quasars themselves (see Figure |2|) . In §H| 
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we show that the excess does not have the same r—K s color distribution as the field population 
described in Appendix |A| (see Figure [TT| ) . In §|] we quantify the strength of the clustering under 
the assumption it is associated with the quasars and look for correlations between clustering 
strength and various quasar properties. In §[| we estimate the amount of luminosity evolution 
bright galaxies have undergone in the K-hand since z=l— 2. In we compare the spectral energy 
distributions of galaxies in several fields with the predictions of various galaxy spectral evolution 
models. We summarize our major results in §||. Some additional details are found in Hall (1998). 



2. Are the Excess Galaxies Associated with the Quasars? 

We consider whether the excess galaxies are plausibly physically associated with the quasars 
by examining the projected radial distribution of galaxies around the quasars and the galaxies' 
distribution in color-magnitude space. An alternate possibility is that they are intervening galaxies 
present because they trace large-scale matter fluctuations which weakly lens the quasars, resulting 
in quasars from radio catalogs with bright flux limits preferentially having galaxy excesses 



around them ("magnification bias"; see Benitez, Martinez-Gonzalez & Martin-Mirones 1997 and 



references therein). We observe 134 galaxies in the magnitude range -fQyx/iiT = 14.5-17 under 
our conservative magnitude scaling, and 150 under our liberal magnitude scaling. The average 
literature counts compiled in §4 of Paper 1 predict 120.3±22.6 galaxies in our survey area, where 
the uncertainty is the RMS scatter among the different surveys in the literature. Thus the excess 
of bright Kukirt<^7 galaxies in our fields is not statistically significant, and magnification bias 
is unlikely to be a significant effect in our sample. 



2.1. Projected Radial Distributions of Galaxies Relative to Quasars 

If the faint excess galaxies in our fields are physically associated with the quasars, they are 
likely to lie preferentially near the quasars on the sky. To test this we examine the projected 
radial distribution of galaxies around the quasars, despite its insensitivity to clusters or groups at 
the quasar redshift but not centered on the quasars. A matched filter ( Postman etal. 1996j ) or 



cell count ( Lidman &: Peterson 1996| ) technique could be used to search for off-center clusters, but 



would be difficult given the small size of our fields. However, the galaxy excess over each entire 
field should still reflect the presence of off-center clusters. 

For this analysis we excluded the fields of Q 0736—063 (uncertain stellar contamination) 
and Q 1508—055 (no r data). We took all galaxies detected at >3cr down to the average 5cr K 
magnitude limit in the 31 other fields, above which essentially no spurious detections are expected 
(see §3.6.7 of Paper 1). We binned the galaxies in 10" annuli centered on the quasar and divided 
by the area imaged within each annulus. No correction was made for loss of objects due to 
crowding or incompleteness or for stars fainter than our star-galaxy classification limits, and we do 
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not count the quasar host galaxy. None of these effects should significantly bias the radial galaxy 
distribution. The results (Figure |2|) show a clear excess of galaxies within 40" of the quasars. The 
data at 9<40" deviates from the level determined from the #>40" data at the 99.999% significance 
level (4.5(7, assuming a gaussian probability distribution). 

We use the reduced chi-squared xt to quantify the deviation from a uniform radial distribution. 
Here v is the number of radial bins used minus the number of parameters determined from the 
data. We use xt both for the deviation from the literature across the entire field and for the 
deviation of the #<40" data from the #>40" data. Since the Xu test requires binning and does not 
distinguish over- and under-densities in individual bins, we also use the Kolmogorov-Smirnov test, 
which requires no binning but is less sensitive to differences in the tails of the distributions (small 
or large radii) than in the middle. 

To investigate the magnitudes of the galaxies producing the observed excess at #<40", we 
take all galaxies with K s <17 and repeat our analysis. The results are consistent with a uniform 
distribution. This is a strong indication that at least some of the faint excess galaxies in these 
fields are associated with the quasars, since both bright (i.e., z<l) and faint galaxies should be 
involved if intervening lensing galaxies caused the excess at #<40" . 

2.1.1. Dependence on Quasar Redshifts 

To search for any dependence of the galaxy excess on the quasar redshifts, we split the quasar 
sample into z<lA and z>lA subsamples. This is a natural division since the two subsamples 
were observed almost entirely in different observing runs, and so any systematic magnitude scale 
offsets will not affect comparisons within them. Also, the z>lA imaging reaches deeper than 
the z<lA imaging in an absolute sense as well as relative to the estimated brightest cluster 
galaxy magnitude at the average subsample z. Repeating our analysis for low- and high-redshift 
subsamples with equal numbers of quasars does not change our results in either a quantitative or 
a significant quantitative sense. Figure ||| shows the ensemble radial profile for the high-redshift 
subsample, which deviates from the uniform level of the #>40" bins at 99.83% (3.15<7) significance 
(from the K-S test). The data at 6»<40" deviates from the 6>>40" prediction at the 99.99% (3.9cr) 
level. The amplitude and spatial profile of the central excess are both consistent with it being 
composed of galaxy clusters at the quasar redshifts. For the z<lA RLQ fields, we plot only objects 
brighter than 0? 1 45 above the average 5<7 K limit, equivalent to 7.5a detections, to ensure uniform 
detection sensitivity at all radii in spite of different exposure times. Figure |] shows a possible 
overall excess within 30", but no excess in the innermost 10". The overall dataset deviates from 
the 6>>40" prediction at only 65% significance (la) and the 6><40" data at only the 97% (2.2a) 
level, due to the large uncertainties (our z<lA field data is shallower than our z>lA field data). 
The #<10" deficit persists even in the complete catalog of 3a detections in either r or K s and is 
not caused by incompleteness due to seeing or inappropriate dithering. This deficit may indicate 
that the z<lA quasars are not centered in any putative host clusters, but excluding the innermost 



- 6 - 



point does not raise the significance of the #<40" excess above 3a. If the #<40" excess is real, the 
6<10" deficit argues against it being due to foreground lensing galaxies. 



2.1.2. Large-Scale Excess 

The excess galaxies in the central ~40" radius region cannot explain the entire excess 
observed in our fields. In Figure || we plot the radial distribution of galaxies to a fixed limit 
of Kjjkirt=^-5 along with the surface density and ±lcx RMS dispersion derived from the 
average published literature counts. Despite the large uncertainties, there are apparently two 
components to the excess galaxies: the central component (#<40") and a large-scale component 
extending to 0~1OO". (At z~1.5 in our H =75, q =0.l cosmology, 40" is ~Q.3h^ Mpc and 100" 
is ~0.75/i7 5 1 Mpc.) Could this large-scale component be spurious? The RMS uncertainties on the 
combined literature counts are in good agreement with galaxy clustering predictions based on 
single surveys ( parlberg et al. 1997 ; Roche, Eales &: Hippelein 1997| ), so the range shown for the 



literature should be accurate. We do not expect many spurious objects above the 5a limits in our 
fields. We have neglected contamination by stars fainter than our star-galaxy classification limits, 
but also incompleteness corrections to our fields (though not to the literature data) . If we include 
both effects, the galaxy surface density is essentially unchanged to Kjjkirt=^0 arid actually 
increases (by <^3%) to Kukirt=20-5- 

A simple check of the need for a large-scale excess can be made: from Figure ||, the central 
excess at #<40" above the surface density defined by the 6>>40" points is only ~10% of the total 
counts, or 0.04 dex, but Figure [j] shows that our excess counts are ~0.1 dex (25%) above the 
literature average at 19<Kjjkirt<^-5. Thus a large-scale excess of amplitude roughly equal to 
the central excess must exist in our fields. Supporting evidence for the reality of the large-scale 



excess can be found in its r—K s color distribution (Figure 14), which is indistinguishable from 
that of the excess at #<40". 



2.1.3. Ubiquity of the Excess Galaxies 

We consider various subsamples of our data to see how common the excess galaxies are. We 
plot the radial distribution of galaxies in our z>lA fields to various limits in Figure ||. Using 
the average surface density from the #>40" data, the #<40" excess has significance 99.7% (3a) 
at K<20.5, and <85% (<lcr) at K<20 or brighter. Using the average literature surface density 
and RMS, the overall excess (central and large-scale) has significance 99.9% (3.3a) at i^<20.5 or 
K<20 and <95% (2a) at i^<19.5 or K<19. Using our conservative magnitude scale effectively 
moves the literature surface density upward slightly. The overall excess then has significance of 
only 96-97% (2-2. 2a) at K<2Q.5 or K<2Q. Thus for our z>\A fields alone, the central excess 
is only significant if the magnitude limit reaches K=20.5, while the large-scale excess is only 
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significant if the magnitude limit is K=20 or fainter, and then only under our liberal magnitude 
scale. Due to shallower magnitude limits and fewer fields, the uncertainties are larger at z<\A 
(Figure ^). The overall excess compared to the average surface density from the literature is only 
significant at the 95% (2a) level at best, and even less under our conservative magnitude scaling. 
We have confidence in the reality of this large-scale excess despite these low formal significances, 
since the limiting factor in determining the significance is the intrinsically large field-to-field RMS 
of K-b&nd galaxy counts. 

The central excess may be produced by only a few of the fields, but the large-scale excess is 
not. If we exclude the 5 (out of 20) z>lA fields with highest x 2 > the central excess drops below 3a 
significance. Thus it may be produced by as few as ~25% of the fields (a very uncertain fraction). 
However, if we take the 15 z>lA fields which reach K=20 and remove the 5 fields with highest x 2 > 
the resulting overall excess is still significant at the 99.7% (3a) level, as measured by both the x 2 
and K-S tests. Thus the large-scale galaxy excess is not produced solely by the same fields which 
may contribute most of the central galaxy excess. At least ~50% of our RLQ fields contribute to 
the large-scale galaxy excess. Both fractions may be higher since removing fields will reduce the 
S/N and significance (to <3a) even when all fields contribute to the large-scale excess. 

Only three of our fields might have a large-scale excess produced by very low redshift galaxy 
associations listed in NED0 within 18': Q 1018+348, 7'A from Abell 982 (no published redshiftO; 
Q 1221+113, 1?7 from the center of the Virgo cluster and 7A from NGC 4352, and Q 2230+114, 
5'.5 from NGC 7305. Excluding them does not significantly reduce the x 2 °f anv fit- Another 



remote possibility is that faint galaxies associated with the supergalactic plane ( De Vaucouleurs 



1975) might affect our counts. Using NED, we found the average absolute supergalactic latitudes 
to be 31+17° for our 0=1-1.4 RLQ subsample, 27+16° for our 2=1.4-2 RLQ subsample, and 
23+14° for published random-field if -band surveys. Thus even if faint galaxy counts do correlate 
with supergalactic latitude, there should be no systematic offset between our counts and those of 
the published literature data due to such an effect. 



2.2. Correlations with Quasar Properties 

Finally, we examine the dependence of the central excess on various quasar properties using 
the x 2 test. Ordering our sample on absolute V magnitude and splitting it in half, we find a 4cr 
#<40" excess for the more luminous RLQs (My=— 27.21+0.61) but only a 2a excess for the fainter 
RLQs (My=— 25.31+0.77). However, these subsamples have average redshifts 2=1.64+0.27 and 
1.32+0.28 respectively. The central excess is more significant in our z>\A subsample than in our 
z<lA subsample. Thus we cannot say whether the primary dependence is on My or redshift. In 
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fact, the apparent redshift dependence may be due to the fact that our z>lA data reaches an 
average of 2™4 below Kbcg {Kbcg is the estimated brightest cluster galaxy magnitude at the 



quasar redshifts), as opposed to 2™1 below Kbcg f° r our z<lA data (see §4.2). This is in turn 
because our 1995 IRIM data (mostly z<lA objects) has higher backgrounds and poorer seeing 
than our 1994 IRIM data (mostly z>1.4 objects). 

Similarly splitting the sample in half at radio power log P ra d=27.5 W/Hz, we find a >4<r 
8<A0" excess for more powerful RLQs (log P ra( i=27.79±0.21) but only 2a for less powerful ones 
(log P ra d=27.14±0.37). The average redshifts of these subsamples are 1.46±0.32 and 1.49±0.32 
respectively, and the average My are -26.20=L1.40 and -26.26±0.95, so this does seem to be a 
dependence on radio power. 

If we compare the 11 flat-radio-spectrum and 19 steep-radio-spectrum objects, we find a 3.6a 
#<40" excess around the steep-spectrum objects but only a 1.65<r excess around the flat-spectrum 
ones. These two subsamples have very well-matched redshift, M a fe s , and P ra d distributions, so 
again the dependence seems to be on radio spectrum. This is supported by the results of splitting 
the sample between sources with and without strong radio lobes (types FRII or T and types C 
or CE, respectively, in Table 1 of Paper 1), subsamples which are 84% and 25% steep-spectrum 
objects respectively. For strong-lobed sources we see a 3Aa 9<40" excess, while for weak-lobed 
sources the excess is only 2. 6a. 

Splitting our sample according to the quasars' associated absorption properties, we find a 
2.8<7 #<40" excess in 13 objects with associated absorption and a 2.4cj excess in 9 objects with no 
associated absorption. These subsamples are slightly mismatched in My, but well matched in z 
and P ra d- Detecting excesses around steep-spectrum RLQs but not ones with associated absorption 
may be surprising in light of the ^>2a tendency for steep-spectrum RLQs to preferentially show 
associated absorption (see the introduction to Paper 1). However, this tendency is not strong 
enough and our subsamples are not large enough to show a significant contradiction. 



2.3. Summary and Discussion: Radial Profiles 

We see an excess of galaxies at #<40" from these 31 RLQs compared to the #>40" background 
level. This central excess is significant at the ~99.995% (~4<r) level and consists of K>17 galaxies. 
It is seen at 3-4ct significance in the z>lA fields, but only at ~2<r significance in the z<lA fields 
which are shallower and fewer in number. The central excess at z>lA may be produced by as 
few as 5 of the 20 fields, but it is not due to one or two extreme outliers. There is an additional 
large-scale galaxy excess extending to #<-^100" detectable in our z>lA fields with 3.3c significance 
at K<20 or _fT<20.5. It may also be present in our z<lA fields, but the larger uncertainties 
there preclude a detection above 2a. We consider various possible errors and conclude that the 
large-scale excess is real at the given significance levels, and that it is not produced solely by 
the same fields which may contribute most of the central galaxy excess. The #<40" excess does 
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not depend on the presence of associated absorption, but seems to be stronger (4a vs. 2a) for 
the more radio-powerful RLQs, and may be stronger (3.6<r vs. 1.65a) for steep-radio-spectrum 
(a r >0.5) RLQs than for fiat-spectrum ones. 

Could the excess galaxies be associated with intervening Mgn absorption-line systems? This 
is difficult to quantify this using radial profiles due to the variation in limiting magnitude between 
fields and the low formal significances of the excesses at all but the faintest magnitudes. However, 



we show in § |4.3| that the 8 z>1.4 RLQs with known intervening absorption and the 8 z>1.4 RLQs 
without it show no differences in their inferred richnesses. Thus there is no evidence that galaxies 
associated with intervening Mgn systems contribute significantly to the observed galaxy excess. 

The large-scale (to 6>~100") excess can plausibly be located at the quasar redshifts. At 
z=l. 4-2.0, a 200" diameter circle corresponds to ~1.5h7 5 Mpc diameter, which is not an 
implausibly large scale, especially if it is connected with the #<40" (0.6h^ Mpc diameter) excess 
which might be e.g. a cluster core or group embedded in a larger overdensity. From the PLE 
model redshift distribution of Roche, Eales & Hippelein (1997), we estimate that ~1.5±0.5% of 
field galaxies to K=20 lie within 5z=0.05 of z=1.70. This redshift bin size is large but similar to 
those in which overdensities have been spectroscopically confirmed at high z ( pickinson 1996b ; 



Steidel et al. 1998a). Our observed galaxy surface density to X=20 is 1.23±0.14 times the 
literature (Figure Thus if all this excess was at the quasar redshifts, it would consitute a 
typical galaxy number overdensity of ~14±10 compared to the model redshift distribution. This 
overdensity is again plausible: virialized clusters are overdense by similar factors on similar scales, 
the RMS galaxy fluctuation on 8/j -1 Mpc scales is os— 1 ( |Lin etal. 1996 ), and superclusters 



have overdensities ~5-40 on ~30/i 1 Mpc scales ( [Small etal. 1998 ). The candidate large-scale 



structures in these fields are not unprecendented; they would be similar to the galaxy overdensities 
of ~10 spectroscopically confirmed by Deltorn etal. (1996; 1997) and Dickinson (1996b) at and/or 
near the redshifts of three z~l radio galaxies on very similar spatial scales. In particular, the 
overdensity near 3C 324 ( Dickinson 1996b| ) is composed of two clumps or sheets of galaxies 



separated by ~7500 km s _1 in their rest frame, evidence for the existence of large-scale structures 
around at least some z>l radio- loud AGN. Similarly, the 8'-separation z=0.6 quasar pair 3C 345 
and Q 1641+3998 may be embedded in a large-scale galaxy structure ( Ellingson fc Yee 1994| ). 



3. Color-Magnitude Diagrams 

Color-magnitude diagrams (hereafter CMDs and denoted as magnitude/color, e.g. K s /r—K s ) 
can provide useful information on the excess galaxies in our fields. Since field galaxies contaminate 
CMDs at all magnitudes and colors, comparison of quasar-field and control-field CMDs must be 
made. We discuss the literature control field datasets we use in Appendix [A|. 

To reduce the uncertainties on the colors, we calculate them using FOCAS isophotal 
magnitudes. We use the same aperture in all filters, namely the isophotal aperture from the 
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coadded r+K s or r+J+K s detection image for the field. This also allows for more accurate 
color limits, since measuring colors through smaller apertures allows very faint objects to be 
detections instead of 3<r upper limits. Comparison of isophotal- and total-aperture colors showed 
no systematics and a scatter consistent with photometric uncertainties. 

Figure || shows the K s /r—K s CMD for stars in 31 RLQ fields. Typically the robust star-galaxy 
separation limit is K c i ass u m =17 .5-18, hence the dropoff in star counts fainter than that. Stars 
at K s >20 were classified from snapshot HST images. Figure || shows the K s /r—K s CMD for 
galaxies in 30 RLQ fields, with upper and lower limits and error bars omitted for clarity. A galaxy 
with a flat spectrum in f v has r—K s =2.06 with our adopted zeropoints (see Appendix A of Paper 
1 and Djorgovski etal. 1995). 

We have statistically corrected the number-magnitude relation at K s >K c i ass n m for unidentified 
faint stars, but how do we statistically account for them in comparing quasar-field and control-field 
CMDs? Stars have a bluer mean r—K s color than galaxies, and thus could bias comparisons 
between fields at different Galactic latitude (and longitude). We could assume that stars uniformly 
populate the color range v — K s — 1—4.5. However, our Kciasslim values are faint enough that stellar 
contamination is only a ~5% effect at fainter magnitudes. We thus make no correction for stellar 
contamination in our consideration of CMDs, but will discuss the effects of it where relevant. 

3.1. One-Dimensional Histograms 

One simple method of comparing CMDs is to split them into different magnitude bins and 
compare the color histograms in each bin. In the following discussion we exclude the fields of 
Q 1508—055 (no r data), Q 2230+114 (nonphotometric Rc data), and Q 0736—063 (uncertain 
stellar contamination). At the 0^5 binning size we use, the ~O7 1 08 difference between our 
conservative and liberal magnitude scales is not significant, so we adopt the latter. In Figure [ll] 
we plot the surface density of galaxies, binned by r—K s color, from all 30 good RLQ fields (solid 
lines) and in the combined opt-IR control fields (dotted lines) discussed in Appendix [A|. Smaller 
histograms represent those galaxies with lower or upper limits to their colors. At Kjjkirt^18, 
there is a clear excess of red galaxies with r—K>5 in the quasar fields. This color and magnitude 
range is where we would expect to see bright cluster ellipticals associated with the quasars. 

To quantify the significance of any difference between the color distributions regardless of 
surface density, taking into account upper and lower limits, we normalize all histograms separately 
to unit sum and use the Peto & Prentice Generalized Wilcoxon Test as implemented in the 
IRAF0/STSDAS task statistics. twosampt. The two r—K distributions are different at only 
67% (la) significance at Kukirt=17-18, but at 99.98% (3.7<r) significance at Kjjkirt= 18-19 



4 The Image Reduction and Analysis Facility (IRAF) is distributed by National Optical Astronomy Observatories, 
operated by the Association of Universities for Research in Astronomy, Inc., under contract to the NSF. 
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and at >99.995% (>4<r) significance in the faintest two bins shown, and at Kjjkirt>21 (not 
shown) for the handful of fields that reach that deep. Thus the excess galaxy population in 
these RLQ fields has a redder r—K color distribution than the field population at >99.995% 
(>4cj) significance. Random field to field variations in r—K color distributions are possible due 
to variations in galaxy population with environment, but should average out over many fields. 
Thus the excess population of predominantly red galaxies must somehow be connected with the 
presence of RLQs in our fields. 

The red galaxy excess persists at Ki/kirt=20-21, in which magnitude range there is a 
significant deficit of blue galaxies. This apparent deficit is difficult to understand, since the 
Kukirt^O.5 faint blue galaxy population is mostly at z<l and should be uncorrelated with 
the presence of z>l quasars in these fields. If foreground galaxy structures cause magnification 
bias in these fields, one might expect slightly fewer blue (late-type) galaxies in the fields due to 
the tendency of red (early-type) galaxies to preferentially inhabit denser environments, but this 
tendency is too weak to explain a complete absence of blue galaxies or why it occurs only in the 
faintest magnitude bin. In addition, the Butcher-Oemler effect ( Butcher fc Oemler 1984] ) is seen 



in clusters to z~l (Rakos & Schombert 1995; Lubin 1996| ), so it is natural to expect that these 



candidate 2=1— 2 clusters should also have a large (~40%) fraction of blue galaxies. 

However, this apparent deficit is also difficult to explain away as spurious. There is adequate 
control field data (~200 galaxies) at Kukirt=2Q-21 from our control fields and the Hubble Deep 
Field IRIM (HDF-IRIM), Moustakas, and Djorgovski datasets. It is possible that conversion of the 
HDF-IRIM and Moustakas V and I magnitudes to r causes part of the offset between quasar- and 
control-field histograms, but the full ~CP5 offset is inconsistent with the scatter in the conversion 
(Eq. g) and the good agreement between control-field datasets obtained with different optical 
filters (Figure |T^). More subtle errors involving different instruments and different magnitude 
measurement techniques cannot be ruled out, but the two largest control field datasets in this 
magnitude range were both obtained on the KPNO 4m with IRIM and catalogued using FOCAS 
total magnitudes and isophotal colors. Using our conservative magnitude scale reduces the deficit 
only slightly. The apparent deficit is of objects with r=22-24 and Kjjkirt=20-21, which is 
well above the detection limits for the r data but within ~0™5-l™0 of our 3a K s limits. The 
Eddington bias ( Eddington 1913| ), the systematic overestimate of faint object fluxes due to the 



increase in number counts with decreasing flux, should only bias our K s magnitudes brightwards 
by ~CPl-(K n 2 at our 5a K s limits, using Eq. 7 of Hogg & Turner (1997). However, this error will 
occur in our control fields as well, and they still show an excess of blue galaxies compared to our 
quasar fields. Any further systematic errors in our magnitude measurements are likely to exist 
only at our faintest limits (i.e. affecting only the K s magnitudes of any "missing" objects) and to 
underestimate the magnitudes of faint objects (e.g. if FOCAS total magnitudes underestimate 
the flux for the lowest S/N objects compared to brighter objects). However, any such errors in K s 
would produce a bias toward bluer r—K s colors, the opposite of what is observed, and would also 
occur in our control fields as well as quasar fields. Lastly, there is no significant deficit of galaxies 
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blue in J-K at £T=19-20 or if=20-21 in our data compared to the HDF-IRIM control field 
data. This suggests that the error occurs either in r or in both J and K s if it is a straightforward 
magnitude measurement error, but we suspect it is instead due to a combination of errors. In any 
case, the excess of red galaxies at Kjjkirt=20-21 is still significant even if we arbitrarily shift the 
quasar- and control-field r—K color histograms so that they match at the blue end. 

In Figures |12| and [13| we plot the surface density of galaxies in each color bin for the z<1.4 
and z>lA subsamples separately. The z<lA subsample is noisier, as usual, and is only useful 
down to Kukirt=20- The red galaxy excess is definitely present at Kjjkirt= 19-20, and perhaps 
at K{jkirt= 18-19 as well. The z<lA and control field distributions are different at the 99.71% 
(3a) significance level at Kukirt=18~19 and at the >99.995% (>4cr) level at Kukirt=W-20. 
The z>lA and control field distributions are different at the 99.99% (3.9a) significance level at 
K UKIRT =18-19, and at the >99.995% (>4a) level at K UKIRT =19-20 and K UKIRT =20-21. 



In Figure 14 we plot the fractional r—K color distributions for galaxies at #<40" (solid 
line) and 8>40" (dotted line) from the quasars. There is no significant difference between the 
histograms. The similarity between the color distribution of the excesses at small and large radii 
suggests that the large-scale galaxy excess is the same population as the 8<40" excess. 

As a check on the reality of the redness of the excess population, we can plot the 



Kjjkirt / J—Kjjkirt CMD for the five fields with J data. Figure 15 shows the K/J—K control 



field CMD (see Appendix 0), and Figure || shows the K/J-K CMD of the five z>\A RLQ 
fields. Comparing the two, there appears to be an excess of galaxies in the RLQ fields with i^;>19 
and J—K^2. In Figure ^ we plot the normalized J—K color histograms in four magnitude bins. 
The J—K distributions are not significantly different at K<19. They are different at the 98.28% 
confidence level (2Aa) at K=19-20 and at the >99.995% (>4a level) at K=20-21, in both cases 
due to the tail at J—K^2 in the RLQ fields. Only the HDF-IRIM dataset has J data for objects 
at K>20. (Note added in proof: the recent data of Bershady, Lowenthal & Koo (1998) has a 
J—K distribution at K s =2Q-2\ consistent with HDF-IRIM.) Thus there is some evidence that the 
excess galaxy population is redder than the field population in J—K as well as r—K. This more 
subtle difference is not surprising since even to K=20 most galaxies are at z<2 QCowie 1996| ) and 
the red envelope of J—K colors for early-type galaxies increases only slowly out to that redshift. 



3.2. Summary: Color-Magnitude Diagrams 

The r—K color distribution of the faint excess galaxy population is significantly redder than 
that of the field population. The J—K color distribution for the five quasar fields with J data 
also shows a red tail not present in the field. There is no significant difference between the color 
distributions of the #<40" and large-scale excess components. The colors and magnitudes of both 
excess populations are thus consistent with a population of predominantly early-type galaxies at 
the quasar redshifts. There is an apparent deficit of blue galaxies in the faintest magnitude bins 
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which is difficult to understand as either a real effect or a single systematic error. However, the 
red galaxy excess is still significant if we arbitrarily adjust the r—K color histograms so as to 
eliminate the apparent deficit. 



4. Estimates of Cluster Richnesses 

We have seen that there is a significant faint galaxy excess in our RLQ fields with a color 
distribution redder than that of the field population. While only spectroscopy can determine the 
redshift distribution of the galaxy excess for certain, its magnitude, color, and spatial distributions 
are consistent with it being composed of clusters or other large structures at the quasar redshifts. 
We now assume that this is the case in order to quantify the richnesses such clusters would have. 

Since we have no information on the dynamical state of these structures, our use of the term 
"clusters" can be somewhat misleading. Spectroscopy is needed to confirm or deny the hypothesis 
that these structures are virialized or virializing. Even then, the evolution of individual galaxies 
and clusters from z>\ to z~0 might make accurate comparison with low-redshift clusters difficult 
without detailed comparison with simulations to identify similar populations of objects at each 



z. We discuss these and similar issues further in §4.3, but for now we assume the strength of the 



galaxy excess in our fields can be compared more or less directly with that of low-redshift clusters. 

Abell (1958) defined a cluster's richness as the number of member galaxies above background 
no fainter than 2™0 below the third brightest cluster galaxy and within 1.5/i _1 Mpc of the center 
center. This criterion is difficult to apply directly at high z due to the large and uncertain 
background correction required over 3/i _1 Mpc diameter and also to the uncertainty in identifying 
the third brightest cluster galaxy. Thus we use the alternate richness measurement iVo.5. 



4.1. The Hill & Lilly Statistic JV . 5 

Hill & Lilly (1991) defined the quantity iVo.5 to be the number of galaxies above background 
located within 0.5 Mpc radius of the quasar and with magnitudes between mscG and rriBCG+3, 
where msec is the brightest cluster galaxy magnitude. This is similar, but not identical, to the 

used by Bahcall (1981). To measure iVo.5 in our data, we need to know mBCG( z )- Also, for 
valid comparison with results at low z, we must account for redshift-dependent changes in the 
bright end of the galaxy luminosity function (between vtlbcg and ttibcg+S)- 



4-1.1. What is rriBCG(z)? 



Hill & Lilly (1991) were studying the environments of radio galaxies at z~0.5, and so 
originally defined mBCc{ z ) to be the magnitude of the radio galaxy in each field. We choose to 
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define mBcc{ z ) as the magnitude given by the K—z relation for high-z powerful radio galaxies 
(HzPRGs), supplemented by the K—z relation of known BCGs at z<\ (see § [4.1.2 ). Assuming 



for the moment that HzPRGs are BCGs (see below for a discussion of this point), this has 
the advantage of having "built-in" evolutionary, cosmological, and k- corrections to the BCG 
magnitudes. In addition, the HzPRG population at z=l— 2 does not show much more scatter in 
the K—z relation than reasonably expected from studies of dozens of BCGs in clusters at z<0.05 



QEales etal. 1993 ; |Thuan fc Puschell 19891 ; postman fc Lauer 1995|) 



One concern with estimating the BCG K—z relation from that of HzPRGs is contamination 
of the if-band light either through scattered or direct AGN light or through AGN-induced star 
formation or line emission. Roche, Eales & Rawling (1997) examined 10 z=l— 1.4 6C radio galaxies 
(RGs), and found they have significantly smaller if-band half-light radii than 3C RGs at similar 
z. Thus the brighter K magnitudes of the 3C galaxies are due, at least in part, to their larger 
sizes rather than an increased contribution by AGN-related flux. Best, Longair & Rottgering 
(1997; 1998) present several other arguments in favor of 3C RGs being BCGs. They claim a 
<^15% AGN-related contribution to the i^-band light, and radial intensity profiles well matched 
by an r 1 / 4 law, both consistent with starlight from old populations dominating at K. In addition, 
the combined 3C RG .fT-band profile shows evidence of excess emission at r>35 kpc, which they 
interpret as cD-type halos. A similar excess of K-band emission at large radii in 3C RGs relative 
to a sample of MG radio galaxies was found by McLeod (1994), who however suggested it may 
be due to very nearby (i.e. interacting) galaxies. There is also some direct evidence that at least 
some 3C RGs at z>l lie in clusters, as expected for BCGs. Dickinson (1997) has spectroscopically 
confirmed a cluster around 3C 324 at z= 1.206 from which X-ray emission has also been detected 



(Smail & Dickinson 1995), and Deltorn etal. (1997) have spectroscopically confirmed a cluster 
around 3CR 184 at z=0.996. Extended cluster-scale X-ray emission has also been detected around 
several 3C RGs by Crawford & Fabian (1996a; 1996b). 

A second issue in determining iVo.5 is whether the galaxy luminosity function (LF) down 
to m#cG+3 has the same redshift evolution as vtlbcg (BCGs are about 2™5 brighter than L* . 
If not, then our values of Nq^ may be biased with respect to low-redshift measurements. Best, 
Longair & Rottgering (1998, their Figure 1) show that at z=\-2 the K—z relation for 3C RGs 
runs slightly brighter than the relation for a nonevolving ellliptical normalized to 3C RGs at 
z<0.6. From nearly complete redshift surveys, Cowie (1996, his Figure 2) shows that at all z<2 
the upper envelope of the field galaxy population's K magnitude distribution tracks the same 
nonevolving ellliptical K—z relation (what Cowie plots is the K—z relation of a nonevolving Sb 
galaxy with Mx=-25.8, but the two relations are very similar at z=\-2). In addition, Cowie etal. 
(1996) estimate the i^-band LF and find, within the uncertainties, an invariant and a to 
z=l and consistency with an invariant M* K at z=l— 1.6. These results show there must be some 
evolution in the BCG and field galaxy populations since neither are ever quite as bright at z=l-2 
as the K—z relation for passively evolving single-burst galaxies with zj=5. However, since both 
populations agree reasonably well with a single K—z relation (that of a nonevolving ellliptical), 
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differential evolution between them is likely to be small at z=l-2, and so our iVo.5 values should 
not be strongly biased compared to ones made at low redshift. 



4-1.2. The K—z Relation for Powerful Radio Galaxies at z=l-2 

To use the K—z relation for HzPRGs to define uibcg i n our quasar fields, we must fit a 
function to that relation. We use the data on HzPRGs presented in Figure 3 of Eales etal. (1993), 
which have been corrected to 34.5 kpc apertures assuming h=0.5, the data of Aragon-Salamanca 
etal. (1993), who give Kqit magnitudes measured within a 50 kpc diameter aperture for 19 BCGs 
in optically selected clusters at z<l, and data for the two spectroscopically-confirmed z>l BCGs 
( Dickinson 1995 ; Stanford et al. 19971 ). Fitting the combined data, we find 



K = (17.12 ± 0.06) + (5.10 ± O.U)log(z) (1) 

with standard deviation 0™37. This relation is consistent with Figure 3 of Eales etal.. We adopt 
Equation [l] as describing the magnitudes of BCGs to z=2. 



4.2. Calculating the Hill & Lilly Statistic N , 5 

We assume Hq=75, qo=0.1, and A=0 to calculate O0.5, the angular distance corresponding 
to 0.5/2.75 Mpc at the quasar z, for each quasar. In each field we extract all galaxies inside and 
outside of #0.5 between Kbcg ( as given by Eq. [l]) and Kbcg+3. However, only six fields have ha 
limits that reach Kbcg+3. (Our z>1.4 data reaches an average of Kbcg+Z™^, while our z<lA 
data reaches an average of Xbcg+2™05.) For each other field, we scale the iVo.5 measurement 
down to K$ a to the expected value down to Kbcg+3 using measurements to both limits in the six 
deepest fields, suitably averaged. This eliminates the bias toward lower iVo.5 observed in the raw 
measurements to K$ a in these fields. The uncertainties on iVo.5 are dominated by the uncertainties 
in the background subtraction. There is no way to reduce the uncertainty on the number of 
galaxies within 0.5 Mpc between Kbcg and .Kbcg +3, but the uncertainty on the background 
subtraction can be reduced by observing larger areas. 

The results (for our liberal magnitude scale) are shown in Figure gj. The plotted error bars 
do not include the uncertainties in the corrections for fields which do not reach Kbcg+S- The 
average iVo.5 value lies above zero since we observe an excess near the quasars. However, the 
uncertainties are large: the average iVo.5 for all 31 good fields is 11.8±12.0, for all 20 z>1.4 fields 
is 11.4±12.8, and for all 11 z<lA fields is 12.5±11.0. Using the relations between iVo.5, B gq , and 
Abell richnesses given in Hill & Lilly (1991), these iVo.5 values correspond to richnesses of Abell 
class 0±1, where by —1 we mean the richness of the general field. In other words, the la upper 
limit on the average near-field richnesses of our quasars is Abell richness class 1, not including 
the uncertainty on the correction to Kbcg+3 where needed. By near-field we mean the galaxy 
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overdensity in the central 0.5h 7 § Mpc radius region compared to the region beyond 0.5h 7 § Mpc 
in our own data. This same Abell richness class Oil result was found for 17 2=1-1.4 5C and 



6C radio galaxies by Roche, Eales & Hippelein (1997). In Figure IE we plot the histogram of 
near-field iVo.5 values for all 31 fields. This plot can be directly compared with Figure 9 of Hill 
& Lilly (1991). Like FR II radio galaxies at z~0.5 (their Figure 9a) and RLQs at z~0.5 (their 
Figure 11), our quasars at 2=1-2 show values extending from ~0 to ~40, with an average greater 
than zero (typical Poisson uncertainties on our points are ±2 bins). 

We now consider the far- field richness, which we define as the galaxy over density in the central 
O.5/175 1 Mpc radius region compared to the literature expectation. We calculate these values in the 
same manner as the near-field richnesses, using the average published literature N(m) and RMS 
in the appropriate magnitude range. The uncertainties are difficult to reduce other than possibly 
through a more accurate determination of the random-field RMS and through observing all fields 
to Kbcg+3- The results (for our liberal magnitude scale) are shown in Figure These plotted 
error bars do include the uncertainties in the corrections for fields which do not reach Kbcg+3. 
The average far-field iVo.5 is larger than the average near-field -/V0.5, as expected. The uncertainties 
are still large: the average far-field iVo.5 for all 31 good fields is 25.2±19.1, for all 20 z>lA fields 
is 28.8±21.4, and for all 11 z<lA fields is 16.4±12.5. This roughly corresponds to Abell richness 
1.5±1.5 for the z>lA subsample, and Abell richness 0.5±1.0 for the z<lA subsample. Thus 
while the average richness levels are higher when compared to the literature expectation, the 
uncertainties are also larger due to the large field-to-field RMS seen in the literature. 

The average "near-field" iVo.5 value for our two control fields is 5.6±8.9, calculated assuming 
2=1.201 and 1.315 for them so that their 5a K limits are equal to Kbcg+3. Their average 
"far-field" A0.5 value is — 13.3±10.3 with the same assumptions. This negative value is expected 
since our control fields' galaxy counts are below the published literature average (Figure 9 of Paper 
1). These control field results show no evidence for any systematic bias in our iVo.5 measurements. 



4.3. Summary and Discussion: iVo.5 

Under the assumption that the excess galaxies are all at the quasar redshifts, and that a 
straightforward comparison of iVo.5 values can be done between these quasars' galaxy excesses at 
2=1-2 and clusters at 2~0, the values of iVo.5 we calculate indicate the near-field excess within 
0.5/w Mpc (~65") around the quasars corresponds to Abell richness class ~0=L1, where by -1 
we denote the richness of the field. This near-field excess is calculated with respect to our own 
data at >0.5/iy 5 1 Mpc. We also measure the excess within 0.5hj^ Mpc compared to the expected 
literature counts; this far-field excess corresponds to Abell richness class ~1.5±1.5. The excess 
across our entire RLQ fields (i.e. to beyond 0.5hj~^ Mpc) compared to the expected literature 
counts is presumably part of the same overdensity as the far-field excess. At some #^100" we 
expect that this overdensity will disappear and our counts will match the literature. 
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Could the excess galaxies be associated with intervening Mgn absorption- line systems? Such 
systems would have to be at z^0.9 to explain the colors of the excess galaxies. From Paper 1, 
intervening absorber information is available for only 4 objects at z<lA, but for 16 objects at 
z>lA. There are 8 z>1.4 RLQs with known intervening absorption (14 z>0.9 Mgn absorbers 
with z=1.33±0.25) and 8 z>lA RLQs without it. These two subsamples have <la differences 
in their near- or far-field N0.5 values despite this ~3.5<r difference in the number of intervening 
absorbers along the line of sight. Although the statistics are admittedly small, there is no evidence 
that intervening Mgn systems contribute significantly to the galaxy excess. This is consistent 



with some published evidence that Mgn and Civ absorbers may tend to avoid clusters ( Bechtold 



fe Ellingson "19921 ; |Morris et al. 1993| ; pirn ;son et al. 1994| ). On the other hand, there is also 



evidence for excesses of galaxies at the redshifts of strong Mgn and/or damped Lya absorbers 
from narrow-band imaging emission-line searches ( Bunker etal. 1995| ; Francis, Woodgate & Danks 



1997] |Bechtold etal. 1998| ; |Mannucci etal. 1998] ; [Teplitz, Malkan & McLean 1998| ). 



One simple interpretation of this result is that RLQs at z=l— 2 are located on a large scale 
(^0.75h^ Mpc) within clusters and/or large scale galaxy structures of Abell richness or greater 
(la lower limit). On a smaller scale (^O.S/i^ 1 Mpc) within those structures, the RLQs can be 
located in environments ranging from the field to clusters up to Abell richness 1 (ler upper limit). 

However, as mentioned in §||, the evolution of individual galaxies and galaxy clusters may 
complicate such comparisons. The numerical simulations of Steinmetz (1997) show that present 
day L* galaxies typically have several progenitors at z~3 spread over a few hundred kpc, and 
Sawicki & Yee (1998) find that the inferred total stellar masses of Lyman-break-selected z^2 
galaxies in the Hubble Deep Field are not large enough for them to be the direct progenitors of 
L* galaxies. If we assume that a present-day L* galaxy has on average two progenitors at z~l-2, 
then our richness measurements will be biased high by a factor of two on average. On the other 
hand, early-type galaxies in clusters or other dense structures may be fully formed at higher 
redshifts than field galaxies, while clusters are probably not fully formed at high z and so the 
galaxy excess in our fields may need to be counted within a larger radius for direct comparison 
with low redshift cluster richness measurements. This illustrates the need for detailed comparison 
of high-redshift observational data with numerical simulations in order to relate such objects to 
their better understood low-redshift counterparts. However, the above considerations will not 
change our basic finding that z= 1-2 RLQs are often embedded in large-scale galaxy structures 
and occasionally also in smaller structures the size of present-day groups or clusters. (See §|| for a 
discussion of these structure's richnesses in a broader context.) 

As large as the uncertainties on A0.5 are, the near-field A0.5 values would have had even larger 
uncertainties if we had not designed our data reduction procedure to make use of the edges of the 
fields which have less than the full exposure time. More accurate determination of RLQ richnesses 
at 0=1-2 will require data down to over wider fields (>4'x4'). Such imaging might best 

be done in J or H if their field-to-field galaxy RMS is smaller than in K. 
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5. The K-band Luminosity Function of Candidate z>l Galaxies 

We have shown that an excess of red galaxies exists in our fields, that the color distribution 
of the excess at #<40" from the quasars is indistinguishable from that at 8>40" , and that the 
color and magnitude distributions of the excess galaxies are broadly consistent with them being 
galaxies at the quasar redshifts. Our data are too noisy for us to fully constrain the luminosity 
function (LF) of the excess population under this assumption, but we can put some constraints on 
M^. First, some terminology: Kbcg is the apparent K magnitude for brightest cluster galaxies 
(BCGs), M^ CG is the BCG absolute K magnitude, M* K is the absolute K magnitude of the knee 
in the luminosity function, and K* is the apparent K magnitude of the knee. 

The usual method used to determine the LF is to assume a cosmology and a prescription 
for ^-corrections, use them to determine absolute magnitudes for the excess galaxies, and fit a 
Schechter function to the resulting absolute magnitude distribution. However, since our fields 
reach only ~2™5 below Kbcg on average, we are not sensitive to (and cannot constrain) the 
faint-end LF slope a. Similarly, since we do not know the volume over which the excess galaxies 
are distributed, we cannot constrain the LF normalization (ft*. 

We assume our usual cosmology of Hq=75 km s~ x Mpc -1 , qo=0.l, and A=0. For our 
no-evolution ^-correction we adopt the average of the -KT-band fc-corrections for the 15 Gyr old 
model E and Sa galaxies of Poggianti (1997). At the relevant rest wavelengths (A>6600 A), these 
/c-corrections change very little over redshifts £=1-2 and ages 11-15 Gyr. We count all galaxies 
across each field brighter than the 5a K limit (K^), correct for incompleteness, subtract the 
expected literature counts, calculate absolute magnitudes assuming the excess galaxies are all at 
the quasar redshifts, sum up the individual fields' excesses, and normalize by the area in each 
absolute magnitude bin. We ignore the faintest bin since it contains data from only one field. 



The results for z<lA and z>lA on our liberal magnitude scale are the solid lines in Figure 21 
and Figure respectively, with dashed lines the ±la Poisson uncertainty ranges. In our 
cosmology our assumed Kbcg (Eq. [l]) corresponds to M|? c ' G =-26.56±0.36 (including intrinsic 
dispersion of CP3) at 2=1.4-2, and -26.28±0.31 at 2=1-1.4. Thus the <2a excess at M K <-27, 
which consists of only a few observed galaxies, is almost certainly not at the quasar redshifts. 

We fit the observed excess N{Mk) with a Schechter function, defined as 

4>{M) = 0.920* exp(-0.92(M - M*)(a + 1) - exp[-0.92(M - M*)]) (2) 



(ISchechter 197§ ). Most previous i'T-band LF determinations find a=-1.0 within their uncertainties, 



and Eq. y reduces to 

<P(M K ) = 0.92(ft*exp(-exp[-0.92(K - K*)}) (3) 

To estimate K*, we simply assume various values of K*, generate (ft(Mx) for each, and shift 
this curve vertically by the weighted difference between it and all bins fainter than Mk=-26.5 
since (ft* is unconstrained. Objects brighter than Mx=-26.5 are very unlikely to be at the quasar 
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redshifts. Bins fainter than M^-=-23.5 in the z=l— 1.4 subsample are excluded due to their large 
uncertainties. The x 2 f° r this shifted curve compared to the data is calculated and and the 
procedure repeated for different K* to find the value which yields the lowest x 2 ■ The result is 
M£=-25.25i^ 10 for our z=l-1.4 fields and M^=-24.85^|^ for our 0=1.4-2 fields. The ±la 
limits are the values of K* —Kbcg which yield ±l<r deviant x 2 values. 



5.1. Discussion 



We compare our results with the LF of Gardner etal. (1997): M^-=-23.90±0.10 at z=0.14. 
Gardner etal. measure magnitudes in 10" diameter apertures ( pardner et al. 1996 ), or 22 kpc 



diameter at z=0.14 for our cosmology. This is equivalent to a 3" diameter aperture at z=l-2, 
smaller than the typical total magnitude aperture for our faint galaxies, but comparable to 
the typical isophotal magnitude aperture. The average offset between these two magnitudes is 
~(K n 3±Cfi n 2 (see Figure 4 of Paper 1). Since our magnitudes are probably systematically brighter 
than those of Gardner etal. by roughly this amount, we include this offset and uncertainty in our 
luminosity evolution estimates. Our results are consistent with mild luminosity evolution, or none: 
-1.05l^ 12 magnitudes in 11 RLQ fields with z=1.13 and -0.65±g;^ in 20 RLQ fields with z=1.67. 
We stress that these estimates are made assuming that all excess galaxies in these fields are at the 
quasar redshifts and with a very simple fit to an arbitrarily normalized Schechter function with 
fixed faint-end slope. Nonetheless, a somewhat more empirical estimate of luminosity evolution 
from <P(K—Kbcg) agrees well with our estimate from <P(Mk). In addition, this is roughly 
consistent with passive evolution of stellar populations formed at high z: Poggianti (1997) give 
passive evolution of -1.05±0.10 and -1.33±0.12 magnitudes to the same redshifts for an equal 
mixture of Sa and E galaxies in an Hq=50, qo=0.225 universe. 

If the excess galaxies were consistently located foreground to the quasars, our estimates of 
their absolute magnitudes would be biased systematically bright. To erase the evolution seen in 

the excess galaxies would need to have 2=0.7 in the £=1.13 RLQ sample, and z=1.2±0.1 in 
the z=1.67 RLQ sample. This is more or less independent of the assumed cosmology. We can rule 
out z~0.7 since galaxies at that redshift would not show such a strong tail to very red colors in 
r—K. We cannot rule out z~1.2 since galaxies can be as red as r— iT~6 at that redshift. Thus our 
detection of luminosity evolution in from <P{Mk) is quite dependent on our assumption that 
the excess galaxies are all at the quasar redshifts. However, if the excess galaxies are at 2^0.9 as 
indicated by their r—K s color distribution, the most likely redshift for them to be at is the quasar 
redshift, since there is no evidence that intervening Mgn systems contribute significantly to the 
galaxy excess (§ |4.3| ) and the only other reason for an excess galaxy population in many fields 
would be weak lensing amplification of these ,2=1-2 RLQs, but large scale structure at z;>0.9 
would not efficiently cause such an effect. 



To summarize, in our adopted cosmology observed i^T-band BCG magnitudes match 
no-evolution predictions to within 0^3. Assuming the excess galaxies are at the quasar redshifts, 



-20- 



we estimate -1.051^ magnitudes of luminosity evolution to 2=1.13 and -0.65I ; 55 to 2=1.67. If 
our assumed cosmology is in error, the amount of passive evolution we infer will also be wrong; a 
(7o=0.5 cosmology would reduce our estimated luminosity evolution by ~0? 1 3. 



5.2. Comparison with Previous Work 



Our result is broadly consistent with previous work. Studies of the iT-band LF through 
various means exist at low redshift ( Mobasher, Sharpies &: Ellis 1993 ; Glazebrook et al. 1995 
Gardner et al. 1997| ), in a z=0.3 cluster ( Barger etal. 1995| ), to z~l (Aragon-Salamanca et al 



1993; Steidel, Dickinson & Persson 1994; Cowie et al. 1996]), and at z>l (Aragon-Salamanca et al 



199% see also [Aragon-Salamanca, Ellis &: O'Brien 1996|) . We adopted the Gardner et al. since 



it is the most accurate value to date and lies between those of the other two low- 2 studies. 

From a K<18 field galaxy redshift survey, Elston (1994) suggested that ~l m luminosity 
evolution to z~l in bright early- type galaxies is required, in agreement with the results of the 
Lebofsky & Eisenhardt (1986) survey of a heterogenous sample of ellipticals and radio galaxies to 
2=1. However, Cowie etal. (1996) find from a K<20 field galaxy redshift survey that within their 
±0^5 uncertainties, is invariant to z=l and is consistent with being invariant at 2=1-1.6, 
where they observe a deficit of bright galaxies (Mr-<-25.1) despite >95% complete spectroscopy, 
which argues against a brighter at z>l. From a Mgn absorption line selected galaxy sample, 
Steidel, Dickinson & Persson (1994) also find no evolution in the -ftT-band LF to 2=1. 

The studies of Barger etal. (1995) and Aragon-Salamanca etal. (1993) found no evidence for 
evolution in (within ±CP 1 3-0? 1 4 uncertainties) in clusters since 2=0.3 and from z=0.37 to 
2=0.9, respectively. (The latter study did find evidence for color evolution of the reddest cluster 
galaxies consistent with passive evolution, however.) Also, Aragon-Salamanca etal. (1994) found 
that the LF of Mobasher, Sharpies & Ellis (1993) was consistent (within uncertainties of ±0^45) 
with that of a sample of candidate z=1.3-2 multiple Civ absorbing galaxies. 

However, there is evidence for surface brightness evolution in both cluster and field early-type 
galaxies from the work of Schade and collaborators ( [Schade etal. 1995 ; Schade etal. 1996ab; 
Schade, Lopcz-Cruz & Barricntos 1997; see Schade 1996| for a review). The amount of evolution at 
a given galaxy size is AMb=-z in the rest-frame S-band to 2=1.2. The simplest interpretation of 
the surface brightness evolution is luminosity evolution at a level consistent with passive evolution 
from formation at high z. There is no evidence for differential evolution between field and cluster 
early-type galaxies to z~0.6, although the dense cluster cores are not well represented in the 
cluster sample. The results on cluster ellipticals are supported to z~0.4 by other spectroscopic 
( van Dokkum Franx 1996 ; Bender, Ziegler fc Bruzual 1996] ) and imaging work (|BarrientosJ 
Schade fc Lopez-Cruz 19*96] ; Pahre, Djorgovski fc de Carvalho 1996 ; Dickinson 19971 ). 



Thus there is no clear consensus on the magnitude of luminosity evolution in galaxies to 
2~1 (cf. Zepf 1997; Hudson etal. 1997; Brinchmann etal. 1998; Lilly etal. 1997). Such evolution 
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should be stronger in the blue, which may help reconcile the r+0±0.20 rest-frame Mb evolution 
of Schade etal. with the <CP5 limit on Mk evolution of Cowie etal.. Given the uncertainties, our 
results are not in conflict with previous work, but the excess population in our RLQ fields does 
suggest a brighter M^ at z>l than at z=0, in contrast to the trends seen by Cowie etal. (1996) 
and Ar agon- Salamanca etal. (1994), but in agreement with the work of Schade and collaborators. 



6. Color Pictures 

Although color pictures themselves cannot really be used for quantitative analysis, they are 
extremely useful for visualization, particularly for comparing the colors of galaxies in different 
fields. The IRAF COLOR. RGBSUN task was used to create Sun color rasterfiles for each field, 
which were converted to color PostScript using xv. The red, green, and blue images were the K s , 
summed r+K s (or J when available), and r images. The images were not convolved to the same 
PSF, which leads to occasional spurious effects. The pictures have maximum intensities set to a 
common surface brightness in mag/pixel 2 , so that a given color and brightness represents the same 
color and magnitude in every field, except for edge effects which arise because the normalized, 
constant-rms images were used to create the pictures. 

Two notable fields are shown here. Even these high-quality plots do not really do justice to the 
color images, which can be better viewed on-line at http://www. astro. utoronto. ca/^hall/ thesis. htm\ . 



In both pictures North is at top, East at left, and the bright, blue quasar at center. Figure 23 



shows the field of Q 0835+580 (z=l. 534). There is a clump of very red (r— K s ^5) galaxies around 



the quasar (which is partially merged with a fainter star.) Figure 24 shows the field of Q 1126+101 
(z=1.516). Of particular note are the faint orange- and red- colored galaxies to the WNW of the 
quasar (bright bluish- white object at center). These galaxies have similar, red r—K s colors, but 
the red-colored galaxies are also quite red in J—K s , which makes them candidate background 
galaxies at z^2.5 (see § |7.3j ). Also of note are the yellow-green objects ~1' ESE and ~2' NNE of 
the quasar. These objects have red r—K s (the ESE object is the reddest object in r—K s in any 
field) but blue J—K s , which makes them candidate extreme late-type stars, although the NNE 
object may be extended. 



7. Spectral Energy Distributions of Candidate z>l Galaxies 

In fields where we have data in filters other than just r and K s , it is worthwhile to compare 
the spectral energy distributions (SEDs) of red objects to models of evolving stellar populations to 
see how strongly the galaxy redshifts and/or stellar populations can be constrained. Our approach 
is a qualitative one since we have at most six points with which to constrain the SEDs and since 
many parameters of old stellar populations are uncertain at the ~30% level using current models 
( Chariot, Worthey fc Brcssan 1996 ; Spinrad et al. 19*9^ ) . 
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7.1. The Candidate Group or Cluster Around Q 0835+580 (3C 205) 

There is a fairly distinct clump of objects within ~20"of Q 0835+580. The 19 catalogued 
objects to our 5a K limit within 20'.'5 (excluding one star at 6>=20'/39) yield a surface density of 
~47 galaxies per arcmin 2 . One object is known to be at z=0.236, and one is unresolved on the 
WFPC2 snapshot of this field and is likely a star. Only one other shows up on the WFPC2 image, 
and it is definitely a galaxy. None of the others are likely to be stars since WFPC2 snapshots can 
detect unresolved objects to very faint limits. The compact spatial distribution of these galaxies, 
and the red color of many of them, makes it very likely that they are at the redshift of the quasar 
(z=1.534) or the intervening Mgn systems (z~1.4368). 

The two reddest objects at 8<2tf!5 have 3a lower limits to their r—K s colors of 7.03 (#339, 
#=11"25) and 6.79 (#347, 6>=10'.'60; this object is in fact an riz dropout with ;>24™9 in each filter, 
and is detected only in JHK S .) We convert their magnitudes to F v and normalize to the same 
flux in K s . These are plotted as the filled points and error bars in Figure for riz J UK s filters. 

To fit these SEDs, we consider various model elliptical spectra at of the appropriate age for 
the assumed z~1.5 of the galaxies, all normalized to the data at K s . In our assumed cosmology 
the universe was 3.53 Gyr old at z=1.534, and is 11.06 Gyr old at z=0, slightly lower than the 
most recent post-Hipparcos globular cluster age estimate of 11.5±1.3 Gyr ( phaboyer etal. 1998| ). 

Exponentially declining star formation histories can adequately fit ellipticals at z=0, but 
cannot fit these galaxies if they are at the quasar redshift of z=1.534. This is shown by the 
dashed line in Figure [25|, which is a 3.4 Gyr old elliptical (E) with an exponentially declining star 
formation rate with r=l Gyr from Poggianti (1997). It fits the data at A>4000/mi in the rest 
frame, but is much too blue at A<4000/iin. We consider models with a shorter star formation burst 
for a better fit. Specifically, we use the Bruzual & Chariot (1996) GISSEL96 spectral synthesis 
code to produce the spectrum of a solar metallicity E viewed 3 Gyr after the start of a 1-Gyr 
burst with a Scalo IMF from 0.1 to 125 Mq, using the theoretical stellar spectra of Kurucz and 
Lejeune, Buser & Cuisinier (1996). However, this model (dotted line) still produces too much flux 
at observed z and possibly i (<4000 A rest), even though it fits well at A>4000 A. (The plotted 
error bars on the fluxes do not include systematic uncertainties, and so we attach little weight to 
the models' not matching the JHK data simultaneously.) As a check on the uncertainties in the 
stellar population modeling, we also use the PEGASE ( [Fioc Rocca-Volmerange 1997|) code to 



simulate a 1-Gyr burst solar metallicity E with Scalo IMF from 0.1 to 120 Mq. Our (qualitative) 
results are unchanged if we use the PEGASE model elliptical, since for this spectral range and 
age it is essentially indistinguishable from the GISSEL model. Different stellar population models 
can often produce very different results (§7), so it is reassuring that our basic conclusions remain 
unchanged whichever of these two models we use. 

To produce a stronger 4000 A break to match the observations, an older or a more metal-rich 
population is needed. Dust extinction will not produce the desired effect because the allowed 
reddening between 0.4-0.9^m rest frame is smaller than the reddening needed at <0.4//m to 
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match the model with the data. Figure 26 shows that a 3 Gyr old 1-Gyr burst GISSEL model E 
reddened by E(B— V)=0.25 (tv~0.7) predicts only ~60% of the observed flux at J. The strength 
and abruptness of the break between z and J (observed) requires a strong 4000 A break between 
those filters, either from age or metallicity. 

A somewhat more speculative possible cause for the strong 4000 A break is a nonstandard 



initial mass function (IMF) with a high low-mass cutoff m£~l-5 Mq (Bithell & Rees 199C; Chariot 



etal. 1993). Depending onmj,, 4000 A break amplitudes ~20-40% larger than those of normal 



ellipticals can occur over a period of ^0.5 Gyr for bursts containing 25-100% of the underlying 



galaxy mass. There is some evidence for high-m^ IMFs in some local starbursts ( [Engelbracht 



1997 ; |Charlot etal. 1993 and references therein). Determining whether a strong 4000 A break is 



caused by a nonstandard IMF rather than high metallicity would require high S/N spectroscopy. 



A 6 Gyr old 1-Gyr burst GISSEL model E (solid line in Figure 25) fits the data reasonably 
well at all observed wavelengths, but this is uncomfortably old for z=1.534. A ~5 Gyr old model 
with burst duration 0.1 Gyr fits similarly well, but 5 Gyr is still quite old for z=1.534. Worthey 
(1994) showed that two populations will appear almost identical if AT/AZ=3/2, where AT and 
AZ are the age and metallicity changes, in percent. In other words, the +100% change in the 
age of the universe at z=1.534 in our cosmology which we require to match the observed SEDs 
of these two objects could instead be produced by approximately a +67% (+0.21 dex) change in 
the metallicity. To test this, we use GISSEL96 to produce spectra for ellipticals formed in 1-Gyr 
bursts with metallicities Z=Z Q ([Fe/H] =+0.0932) and Z=2.5Z Q ([Fe/H]=+0.5595). The results 



are shown in Figure 27. The high-metallicity E (solid line) is as good a fit to the data as the 6 Gyr 
E model in Figure 25. Thus the SEDs of these objects are consistent with ~3 Gyr old metal-rich 
ellipticals, even though the metallicities of z=0 ellipticals averaged over the entire galaxy like our 
SEDs are likely to be roughly solar ( |Spinrad et al. 1997 and references therein). 



However, these two objects are redder than LBDS 53W091 (z=1.552), a mJy radio source 
with r— iT~6.15 ( Dunlop et al. 1996| ) whose ultraviolet spectrum is remarkably similar to M32. 



M32 is believed to contain an intermediate-age stellar population ~4-5 Gyr old in addition to the 
old stellar population usually present in ellipticals ( [Spinrad et al. 1997 , and references therein). 



Spinrad etal. (1997) claim an age of ^3.5 Gyr for LBDS 53W091, and Dunlop (1996) quote 
^4.5 Gyr for a similar galaxy (LBDS 53W069 at z= 1.432). On the other hand, Bruzual & Magris 
(1997) claim that LBDS 53W091 can be modelled as an elliptical only 1-2 Gyr old if the full 
rest-frame 0.2— l//m SED is taken into account, The discrepancy is partially due to the evolutionary 
spectral synthesis models they use ( [Bruzual fc Chariot 199 j ), which produce red rest-frame R—K 
colors as fast or faster than any other such models ( [Spinrad et al. 1997 , Figure 16). However, it 



is also partially due to the rest-frame optical-IR colors of LBDS 53W091 being bluer than those 
of M32 despite the similarities in their rest-frame UV spectra. (Heap etal. (1998) also claim an 
age of ~2 Gyr for LBDS 53W091, on the independent basis of comparison with a Space Telescope 
Imaging Spectrograph ultraviolet echelle spectrogram of an F8V star.) 
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In Figure |28| we compare the SEDs of our two r—K s ^7 objects to the spectrum of M32 
( Bica etal. 1996| ), normalized at K s . Within the systematic uncertainties, the M32 spectrum is 



a good fit. This is consistent with the good fit of the 6 Gyr model in Figure 25. The Bica etal. 
M32 agrees well with the 0.4-0.65/im spectrum of Hardy & Delisle (1996), but is bluer than the 
"M32(b)" SED of Magris & Bruzual (1993). Even adopting the latter spectrum, however, the 
main objection of Bruzual & Magris (1997) to the old age claimed for LBDS 53W091 at z~1.5 
is substantially weaker for these two galaxies at similar z, and possibly for a third object with 
r—K s ^6.5 in the field (#618, 9=99"). Deep optical spectra of these objects have the potential 
of placing constraints on the cosmological model as strong as or stronger than those claimed for 
LBDS 53W091 by Spinrad etal. (1997). 

The seven objects next reddest in r—K s at 9<2G/.'5 have r—K s =4.95-5.80. The remaining 
nine at 0<2Cf.'5 — excluding the galaxy of known z=0.236 — have r— X s =2.68-4.49, consistent 
with any z^0.5, and we do not consider them further. We normalize these objects to the same flux 
in iT s and compute the weighted average F v in each filter. We plot the average and RMS scatter in 
Figure ^ and compare to model spectra. The 3.4 Gyr old Poggianti E model produces too much 
flux in observed riz. A 2 or 3 Gyr old GISSEL model is a decent fit: the 3 Gyr old model (dotted 
line) fits best in z, and the 2 Gyr old model (solid line) fits best in i, but both predict too little 
flux in r. This can be interpreted to mean that, on average, star formation dropped off in these 
objects more quickly than an exponential but less quickly than a sudden cutoff of a burst. The 
inferred ages of these very red objects within 20" of 3C 205, and others at larger distances from the 
quasar, do not strongly constrain cosmological models. However, if they truly are younger than 
the extremely red objects (EROs) with r—K s ^6.5 discussed earlier and yet are also early-type 
galaxies at the same redshift, it would mean that z~1.5 is near enough to the formation epoch of 
these galaxies for the dispersion in the rest frame U—V colors of early-type galaxies to rise above 
the small values seen at z<^l ( Stanford, Eisenhardt Dickinson 1998j ). 



7.2. The Southern Radio Hotspot of 3C 205 

Lonsdale & Barthel (1986; 1988; 1998) have studied the radio morphology of this quasar at 
VLBI resolution. The southern hotspot is extremely compact, curved, and polarized, and is is 
located 2"3 W and 8'.'2 S of the quasar. One of the extremely red galaxies discussed above (#347) 
is located 3'/80 W and 9 / /89 S of the quasar (0=lO / . / 6O), consistent with the hotspot being produced 
by the collison of the radio jet with part of this galaxy. (There may be a small (<3°) rotation of 
our image relative to true north.) This "hotspot galaxy" has -ftT s =19.53±0.07, r— K s >6.79 (3a 
lower limit), z-J>3.89 (3a), J-if s =1.87±0.06, and H-K s =0.85±0.18. As we saw in §[?U], this 
object's SED is red enough to match that of a metal-rich and/or very old galaxy, but its observed 
>l/mi slope does not allow much dust extinction. Thus the jet is probably being deflected by a 
high-pressure X-ray halo around the galaxy rather than dusty gas within it. 

Could the hotspot object be synchrotron emission from the hotspot itself? The SED of such 
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emission is expected to have f v oc v~ l ( |Ridgway fc Stockton 1997 ), or r—K s =3.33, z—J=lA9, 
J—K s =1.61, and H— K s =0.81. The near-IR colors of the hotspot object are consistent with the 
expected synchrotron emission, but such emission cannot explain the break at ~l/im, observed, 
which is most naturally explained as a 4000 A break at z~1.5, consistent with the quasar redshift. 



7.3. Candidate Associated and Background Galaxies in the Field of Q 1126+101 

A very intriguing feature of the field of Q 1126+101 is that many of the r— K S >A galaxies 
have very red J—K s colors as well (cf. the reddest objects in Figure 24). In this section we show 



that the SEDs of these objects are most consistent with either galaxies at £^2.5 {background 
to the quasar) whose 4000 A breaks lie between J and K s , or with very dusty galaxies at the 
quasar redshift of z=1.516. Many of these objects have r—K s and J—K s colors similar to the 
spectroscopically confirmed z=2.38 galaxies of Francis, Woodgate & Danks (1997), and several are 
actually redder in J—K than the only ERO with a spectroscopic redshift, HR10 (|Hu fc Ridgway 



1994 ), which requires dust reddening in any reasonable cosmology to fit its SED at its known 



z- 



1.44 QGraham fc Dey 1996; ; |Cimatti etal. 1997| ). 



We attempt only qualitative SED modelling, first trying to fit each object with elliptical 
spectra of reasonable age for the assumed z to minimize the required dust extinction. To model 
the effects of dust, we used the standard Milky Way (MW), LMC, and SMC extinction laws and 
the Calzetti (1997) law. The Calzetti formula was derived for active star- formation regions and 
empirically incorporates the "selective attenuation" effects of dust, including extinction, scattering, 
and geometrical dust distribution effects. In any case, the relative extinction from ~2500 to 
9000 A (observed r to K at #=1.5) is similar for all four of these curves, as seen in Figure 30. 

In Figure [}(] we plot the two objects with r—K s >5 which are reddest in J—K s (#424 and 
#315). A 2.2 Gyr old r=l Gyr Poggianti E at z=2.5 with E(B-V)=0.35 (iy~l) provides a good 
match to the SEDs of these objects. The choice of reddening law does not strongly affect the 
fit except at r and possibly at z. The break between J and H is suggestive of z~2.5 for these 
two objects, but the lack of detected flux except at r and K s makes for a weak constraint on the 
redshift. However, a lower redshift would require a larger ry since longer rest-frame wavelengths 
would be observed, and a younger galaxy or one with a more extended star-formation history 
would require more dust. The required reddening can be reduced to E(B— 1/)~0.25 by fitting with 
a ~2-Gyr old 1-Gyr-burst and a young (^0.1 Gyr) population contributing substantial light only 
at <l/iin (observed). Similar E(B—V) values were found by Sawicki & Yee (1998) for confirmed 
z>2 galaxies in the Hubble Deep Field (HDF) through detailed SED fits. 

Figure [H] shows that the same E(B— V)=0.35 Poggianti E fits the J and K s fluxes of the 
two objects with r—K s <5 which are reddest in J—K s (#18 and #237). However, these objects 
definitely require an additional young population to match the blue r—J colors. 



The lack of usefully deep H data for the above objects prevents firm redshift estimates. Our 
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assumed z=2.5 is the lowest that places the 4000 A break redward of J, but while strong breaks 
in the objects' spectra are consistent with the observed SEDs, they are not necessarily required. 
Figure |3^ shows the only two objects with J—K s >2.5 with firm detections in H which might help 
better constrain their redshifts (#425 and #381). These objects are more consistent with dusty 
galaxies at the quasar redshift z=1.516 than objects at z~2.5 due to the strong spectral break 
between z and J (compare to the objects in Figures and 31 which cannot have such strong 
breaks). Such objects are unusual; there are only two other J— K s >2.5 galaxies which definitely 
have similarly strong breaks. The red J—K s and blue r—z colors can be fit only with composite 
stellar populations, e.g. a dusty 2 Gyr old population with E(B— V)=0.75 (ry~2) and a young 
(100 Myr) dust-free population with ~1% of the total mass. The only way to avoid dust would be 
with metallicities Z>5Zq and/or ages ^5 Gyr, plus a young component providing flux in r. 



7.4. Candidate Dusty Associated and/or Background Galaxies in Other Fields 

The presence of candidate dusty galaxies around RLQs at z~1.5 is relevant to the evolution 
of early-type galaxies in clusters and groups, and the surface density of candidate 2^2.5 galaxies 
is interesting in its own right, so we briefly examine J—K s >2.5 objects in the other fields with 
deep J data (Q 0835+580, Q 0952+179 and Q 1258+404). 

The surface densities of objects with J—K s >2.5 and K<21 are 4.31 arcmin -2 in the 
Q 1126+101 field and 1.90 arcmin -2 in the other three fields. If we also require r<25.5 the surface 
densities are 2.80 and 1.29 arcmin -2 respectively, to ±25% accuracy. Whatever these galaxies are, 
the Q 1126+101 field is clearly unusual. 

It is worthwhile to compare the surface density of these candidate z^2.5 "J-dropouts" 
(objects with J-K>2.b) with that of "[/-dropouts" at z~3.0±0.5. From Steidel etal. (1998b), 
to i?AB=25 (or r=25.5) there is ~1 [/-dropout/arcmin 2 . The four [/-dropouts with published 
K data have K=2l. 3-22.1 and ^5-^=2.7+0.5 ( [Steidel et al. 1996|) . Thus, at the same r 
magnitude the surface densities of J- and U -dropouts are approximately equal, but our J-dropouts 
are about ~2 m brighter in K and thus ~2 m redder in r—K than [/-dropouts. This implies that 
many J-dropouts should be found among [/-dropouts studied in the IR. Only four with J—K>2.5 
have been seen among the four objects mentioned above and the 17 [/300-dropouts in the HDF 
studied by Sawicki & Yee (1998). These results could be reconciled if faint z~2.5 galaxies are on 
average redder than faint z~3 galaxies, but if anything the opposite is seen by Sawicki & Yee 
(1998). To reduce the discrepancy between the inferred surface densities of J-dropouts at 2^2.5 
and [/-dropouts at z~3, we are led to the conclusion that a substantial fraction of our J-dropouts 
are not at z~2.5 but are dusty galaxies at the quasar redshifts. Even then, the relative surface 
densities are such that we expect up to 50% of [/-dropouts will be J-dropouts, larger than the 
20+10% seen to date among the 21 [/-dropouts with published near-IR data. 

This conclusion is consistent with the surface densities of J-dropouts in literature control 
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fields with J data (Appendix [AJ). In 43.14 arcmin 2 to K=21 the surface density of J—K s >2.5 
galaxies is 0.46 arcmin -2 , about half that in our fields (excluding Q 1126+101). This suggests 
that roughly half our J-dropouts are in fact associated with the quasars. 

The SEDs of our J-dropouts are also consistent with about half being dusty galaxies at the 
quasar redshifts z q . Only Q 0835+580 and Q 1126+101 have z- and J-band data which is needed 
to identify the 4000A break at z ? ~1.5. If we select K s <21 objects with J— K s >2.5 and r—K s >5 
(indicating a possible break between r and K s ), 7 of 22 such objects in the Q 1126+101 field and 
7 of 10 in the Q 0835+580 field (roughly half in all) have z—J colors or lower limits consistent 
with dusty galaxies at z q . 



7.5. Summary: Spectral Energy Distributions 

There is a distinct overdensity of galaxies around Q 0835+580 (3C 205). The extremely red 
spectral energy distributions (SEDs) of the two reddest objects within ~20" (and a third object 
99" away) could be due to metallicities ^+0.2 dex above solar or to stellar populations as old 
as M32 (4-5 Gyr) at z~1.5. Dust extinction cannot explain the entire SEDs due to the strong 
break between z and J (observed), presumed to be the 4000 A break. The SEDs of the seven 
next reddest objects in r—K s within ~20" of 3C 205, and of many other galaxies across that field, 
are consistent with 2-3 Gyr old ellipticals formed in a 1-Gyr burst. If these galaxies and the 
three even redder ones are all early-type galaxies at the quasar redshift, then we may be seeing a 
dispersion in the colors, metallicities, or ages of early-type cluster galaxies or their progenitors. 

In the field of Q 1126+101, there are three dozen objects with J—K s >2.5, a distinct excess 
compared to other fields, which themselves show an excess of such "J-dropouts" compared 
with literature control fields. Several lines of argument suggest about half of these galaxies are 
most likely old and dusty (~2 Gyr, E{B— F)~0.75) galaxies at the quasar redshifts with young 
(<0.1 Gyr) components comprising ~1% of their masses. The remaining J-dropouts are consistent 
with being at z^2.5, which places the 4000 A break between J and K s . 

The possible presence of relatively old, dusty galaxies at the quasar redshifts coupled with the 
possible dispersion in early-type galaxy ages mentioned earlier suggests there may be considerable 
dispersion in the properties of cluster ellipticals at z~1.5, in agreement with hierarchical clustering 



models (Kauffmann 1996). However, these candidate relatively old, dusty galaxies could also be 



younger systems with higher E(B—V). 

8. Conclusions 

The major results of our study so far are: 
• We find a significant excess of -f^;>19 galaxies in the fields of 31 z=\-2 RLQs, on two spatial 
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scales (§2.1). One component is at 9<A0" from the quasars and is significant compared to the 
galaxy surface density at #>40" in the same fields. The other component appears roughly constant 
across the fields to #~100" from the quasars and is significant compared to the galaxy surface 
density seen in random-field surveys in the literature. The #<40" component may be produced by 
as few as ~25% of the z>1.4 fields but the large-scale component is present in ^50% of them. 

• The r—K color distributions of the excess galaxy populations are indistinguishable from each 
other and are significantly redder than the color distribution of the field population, consistent 
with the excess galaxies being predominantly at z>l (§ |3.1| ). However, there is a deficit of blue 
(r— K<^3.5) galaxies at K=2Q-2\ which is difficult to understand as either real or spurious. 

• The magnitudes and colors of the excess galaxies are consistent with a population of mostly 
early-type galaxies at the quasar redshifts, such as would be found in quasar host clusters or 
groups. There is no evidence that they are associated with intervening Mgn absorbers (§ [4.3| ). 

• The average excess within 0.5/ifg 1 Mpc (~65") of the quasars corresponds to Abell richness 
class ~0±1 compared to the galaxy surface density at >0.5hj^ Mpc from the quasars, where -1 
denotes the richness of the field, and to Abell richness class ~1.5±1.5 compared to the galaxy 
surface density from the literature (§||). This assumes that the excess galaxies are all at the 
quasar redshifts. This suggests that on a large scale (^0.75h^ Mpc) RLQs at £=1-2 are located 
within clusters and/or large scale galaxy structures of Abell richness ~1.5. On a smaller scale 
(^O.S/iyg 1 Mpc) within those structures, RLQs can be located in unremarkable "field" environments 
or in groups or clusters of Abell richness ~0. One uncertainty not taken into account when 
calculating richness measurements is the correspondence between individual galaxies at z>l and 
z~0. A high merger rate at z<l~2 would mean that individual galaxies at z^O were composed of 
several progenitors at z>l-2, which could bias our richness measurements high. Spectroscopy and 
detailed comparison with numerical simulations are needed to resolve this uncertainty. 

• The only significant correlations observed between environments and quasar properties is that 
the galaxy excess at 8<40" seems to be stronger for the more radio-powerful and steeper-radio- 



spectrum RLQs (§2.2). The 8<A0" excess does not seem to depend on the presence of associated 
absorption. These dependences are based on only 31 RLQs (or fewer), but they at least illustrate 
trends which could be verified with larger datasets. 

• By assuming the excess galaxies are at the quasar redshifts and fitting their i^-band luminosity 
function, we find -0.651q55 magnitudes of luminosity evolution in MJ^ to z=1.67 (§[|). This is in 
contrast to the trends seen at z>l by Cowie etal. (1996) and Aragon-Salamanca etal. (1994), but 
plausibly in agreement with the work of Schade and collaborators ( (Schade 1996| ). 

• For four fields with data in at least r,JK s , we find that the SEDs of most of the excess galaxies 
are consistent with them being 2-3 Gyr old early-type galaxies at the quasar redshifts of z~1.5, 
but that there are galaxies whose SEDs cannot be fit by such simple models (§(7|). At least three 
objects in these four fields have SEDs consistent with being 4-5 Gyr old at z~1.5, and at least a 
dozen others are consistent with ~2 Gyr old but dust-reddened galaxies at the quasar redshifts 



-29- 



(§ |7.4|) . Taken together, these potentially different galaxy types suggest that there is considerable 
dispersion in the properties of early-type cluster galaxies at z~1.5, in agreement with hierarchical 
clustering models ( Kauffmann 1996| ). However, these galaxies could also be younger systems with 



higher E(B—V). Spectroscopic followup will be needed to confirm this suggestion. In particular, 
age determinations from deep spectra of the candidate 4-5 Gyr old galaxies offer the possibility of 
constraining the cosmological model by requiring a relatively old universe at large lookback times. 

• There are also several dozen galaxies in the four fields with good J data (particularly in the 
Q 1126+101 field) whose SEDs are best explained either as background galaxies at z^2.5 or very 
dusty galaxies at z=z q (§f7j]). Some of the z^>2.5 candidates seem to be dusty, to have composite 
stellar populations, or both, and some may be already ^2 Gyr old at z^2.5, again offering a 
possibile constraint on cosmological models. 

Without spectroscopic confirmation it is premature to draw firm conclusions about the 
implications of the existence of the excess galaxy population in these z=l~2 RLQ fields. 
Nonetheless, if we assume that the excess galaxies are predominantly at the quasar redshifts, we 
can draw some interesting conclusions which illustrate the value of spectroscopic followup. 

The existence of galaxy excesses on two spatial scales, consistent with Abell richness ~0±1 
clusters embedded in Abell richness ~1.5±1.5 structures if the excess galaxies are at the quasar 
redshifts, might lend support to hierarchical clustering models independent of Vl. This is because 
in such models clusters are built up by accretion of smaller clusters and groups over time. Galaxy 
structures embedded in larger structures such as we observe could occur during some phases of 
cluster formation in hierarchical clustering models. 

This hierarchical clustering interpretation of the galaxy excesses on two spatial scales may 
tie in with RLQ host cluster observations at z<0.7 (see §|l] and the introduction to Paper 1). 
At such redshifts the clusters appear to be younger and less virialized than optically-selected 
clusters ( Ellingson, Green &; Yee 1991| ) and are of Abell richness 0-1, with a few of richness 2. 



At z=l-2 we also find that RLQs can be found in overdensities of Abell richness 0-1, but that 
these "near-field" overdensities are embedded in larger structures of Abell richness ^1. It is an 
open question whether the most important factor affecting the evolution of RLQs in clusters 
is the absolute density of the environment, the density within 0.5 Mpc, or even conditions in 
the host galaxy and any galaxies interacting with it. However, despite the fact that an Abell 
richness (sub)cluster observed at 2=1.5 will evolve into a very different structure at z=0 than 
an Abell richness cluster observed at z=0.5 will, we observe RLQs in Abell richness ~0 (within 
O.5/175 1 Mpc) clusters at both z~0.5 and z~1.5. This may imply that it is the properties of a 
RLQ's environs within 0.5 Mpc which determine whether the environs are hospitable to RLQ 
formation and fueling; whether or not the Abell richness ~0 cluster is embedded in a larger 
structure or not may be immaterial. Alternatively, if clusters embedded in different large-scale 
environments show important differences in their formation histories, then the large-scale density 
of the RLQ's environs would be what determines which clusters host RLQs. If this is the case then 
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the important result is that quasar host clusters at z~1.5 are richer on large scales than at z~0.5. 

If we assume instead that the interpretation of Ellingson, Green & Yee (1991) is correct, 
i.e. that galaxy interactions are largely responsible for creating and fueling RLQs and that only 
dynamically young clusters are hospitable to RLQ, then one interpretation of our z=l-2 results 
would be that we are seeing dynamically young subclusters which will later virialize and/or merge 
with other galaxies or subclusters in the observed large-scale galaxy structures to form clusters 
as we know them today. The fact that RLQ host clusters at z=l-2 are not markedly richer on 
O.bhj^ Mpc scales than those at z<0.7 is again consistent with hierarchical clustering models. In 
such models richer clusters do not undergo monolithic collapse at higher z, but instead form from 
mergers of subclusters at higher z. 

The above speculations should be treated as such. Spectroscopic investigation of these fields 
and detailed comparison to numerical simulations is needed before we can begin to formulate 
a true model of the evolution of quasar environments to z=2. Nonetheless, if and when such 
confirmation is obtained, we feel that the above discussion outlines a reasonable framework for 
future thinking on this subject. 

Among the followup research projects suggested by this work, two will be very important 
in confirming (or denying) and extending our results. First, larger areas (>4'x4') around our 
sample need to be imaged to confirm the large-scale (0~1OO") galaxy excess around z= 1-2 RLQs 
and to determine the true angular extent of this excess. Second, multislit optical and near-IR 
spectroscopy is needed to verify the existence of overdensities at the quasar redshifts, to determine 
velocity dispersions of any large-scale structures found, and to discriminate age, metallicity, and 
dust effects in the brightest candidate z>l galaxies. Over the next decade IR spectrographs will 
begin to appear on 8-m class telescopes (e.g., the MMT 6.5m, LBT, and Keck) and perhaps the 
Next Generation Space Telescope. These instruments will enable spectroscopic identification of 
very red and faint galaxies and studies of their ages, metallicities, and dust reddenings. Ultimately, 
they may enable the fundamental plane of elliptical galaxies and its evolution to be studied to 
2=1-2 and beyond using galaxies from this and other work. 
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A. Literature Control-Field Data 

Since our own control fields cover far less area than our RLQ fields, it is advantageous to have 
as much additional control field data as possible, preferably with data in r and K s . Of the datasets 
available in the literature, however, few utilize those specific filters. For purposes of comparing 
color-magnitude diagrams, it is acceptable to use data in various optical filters to estimate r 
magnitudes, since the uncertainties in predicting r—K s from I—K (e.g.) for faint galaxies are 
smaller than the binning size of CP5 we will typically use in our comparisons. To estimate galaxy 
colors in one optical filter set from another, we use the galaxy colors for morphological types 
E through Im which have been computed by Frei & Gunn (1994; FG94) from z=0-0.6 using 
the SEDs of Coleman, Wu & Weedman (1980), and by Fukugita, Shimasaku, & Ichikawa (1995; 
FSI95) from z=0-0.8 using the spectral atlas of Kennicutt (1992). Specifically, we use Tables 
2-6 of FG94 and Tables 3-9 of FSI95 where both references give consistent first-order color-color 
conversion results (after accounting for different zeropoints in the two references, cf. Appendix A 
of Paper 1) and FSI95 alone where they do not. 

At K s <20, for r—K s comparisons we use primarily our two control fields (18.78 arcmin 2 , Rc 
and K s ), the McLeod etal. (1995) Her-1 field (10.56 arcmin 2 , Tyson R T and K s ), the HDS (|Cowie| 
et al. 1994| ) SSA13, 17, and 22 fields (5.33 arcmin 2 , U'BVI C and K), and preliminary data (25 
arcmin 2 , R c and K s ) from the SA57 field of the 100 arcmin 2 BRIzJK KPNO 4m+IRIM survey 
of Elston, Eisenhardt, k Stanford (1998; EES98), kindly provided by R. Elston. At K s >20, we use 
the Djorgovski etal. (1995) Hercules field (0.44 arcmin 2 , r and K s ), the Moustakas etal. (1997) 
Fields I and II (2.88 arcmin 2 , VIK), and the preliminary KPNO 4m+IRIM catalogs of the Hubble 
Deep Field (HDF-IRIM; 7.40 arcmin 2 , F606W, F8UW, JHK), kindly provided by M. Dickinson 
(cf. pickinson et al. 1998 ). Conservative uncertainties of iO? 1 ! were assigned to magnitudes, 



and ±0™15 to colors, where none were given. We only removed objects classified stellar by the 
various authors, who used varying criteria including morphologies and colors (cf. Saracco etal. 
1997). Stellar contamination is not a large effect at faint magnitudes, but should be kept in mind 
nonetheless, particularly for the HDF-IRIM and EES98 datasets which still include stars. 

We estimated the r and i magnitudes of the McLeod objects from their Rt and It magnitudes 



( Gullixson etal. 1995 ) using relations derived from FSI95 and our adopted zeropoints (see 



Appendix A of Paper 1 and Appendix |A| of this paper): 

R T + 0.367 - 0.048 x (R T - I T ) (Al) 



r 

i = I T + 0.710 + 0.045 x (R T - I T ) (A2) 



We estimated the r magnitudes of the EES98 objects from their Rc magnitudes using a relation 
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derived from both FG94 and FSI95: 

r = R c + 0.322 (A3) 

This is the same relation used for our own control fields and the Q 2230+114 field, and the color 
term is negligibly small. We estimated the r magnitudes of the Moustakas and HDS objects from 
their V and Iq magnitudes using a relation derived from both FG94 and FSI95: 

r = I C + 0.447 + 0.452 x (V - I c ). (A4) 

The lcr scatter around this fit is the largest of all our fits, but is still +0? 1 ! or less at z<^l. 
(The published V and I magnitudes for Field II of Moustakas et al. are in error and should be 
adjusted by —VP 6 and +VP6 respectively.) For the Dickinson etal. (1998) data, we convert the 
HDF F606W and F814W AB magnitudes to Vega-based magnitudes which we dub Vqoq an d 
IgH using V 606 =F606W -0.116 and I 8 i 4 =F814W-0.439 (|Williams et al. 1996|) . We then convert 
to V and Ic using F=V606+0-37(V606 — Isia) and Ic=I%u— 0.10(V606 — hi4:)i which are first-order 
approximations derived from Tables 10 and 7 of Holtzman etal. (1995). Finally, we estimate r 
from V and Ic as for the Moustakas data above. 

Control field J band data were taken directly from the McLeod etal. (1995) SA57SO field, 
the EES98 SA57 field, and the HDF-IRIM data ( pickinson et al. 1998; ). 



We use several different combinations of these datasets in our work. We exclude the Soifer 
etal. (1994) data since those fields were targeted around high redshift objects and are thus 
not the random fields we desire as control fields. By "published control fields" we refer to all 
datasets mentioned above, except the HDF-IRIM and EES98 datasets, plus the bright-end data 
of Glazebrook etal. (1994) and the HMWS and HMWS QGardner 1995a| ; pardner 1995b| ). By 



"opt-IR control fields" we refer to all datasets from which we generated r—K s colors, including our 
own. This includes only parts of some datasets: all but the HDS SSA4 field of Cowie etal. (1994), 
just the Hercules field of Djorgovski etal. (1995), just the Her-1 field of McLeod etal. (1995), and 
just the SA57 field of EES98. 

To check that our conversions of the various literature colors to r—K colors are accurate, we 



compare r—K distributions in Figure IC. The dotted lines are all opt-IR control fields except our 
own, and the solid lines are our own control fields and those of Djorgovski etal. (1995), which are 
the only fields with data originally taken in r (or R) and K s . The uncertainties are large, but 
there are no significant differences between the actual and converted r—K distributions at any 
K=\7-2\, i.e. between our control field data and that from the literature, given the deeper limits 
of many of the opt-IR control fields. 
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Fig. 1. — The conservative Kukirt N(m) relations for our 1994 and 1995 KPNO 4m run data 
are plotted as dotted and dashed lines respectively. The area-weighted average of our conservative- 
magnitude control fields and all published random-field imaging surveys (corrected to Kukirt) is 
plotted as the solid line. RMS uncertainties are plotted for the N(m) values and formal uncertainties 
for the magnitude bin centers. 
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Fig. 2. — Radial distribution of galaxies relative to the quasars. All galaxies detected at >3cr 
down to the average 5cr magnitude limit in each field are included. Error bars are calculated for the 
number of objects in each bin as per Gehrels (1986). The large uncertainties are from small number 
statistics: at small radii the bin area is small, and at large radii only part of the full annulus is 
imaged, so the number of galaxies detected is small in both cases. 
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Fig. 3. — Radial distribution of galaxies in the z>1.4 fields. All galaxies detected at >3a down to 
the average 5a magnitude limit in each field are included. 
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Fig. 4. — Radial distribution of galaxies in the z<1.4 fields. All galaxies detected at >3<r down to 
the average 7.5a magnitude limit in each field are included. 
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Fig. 5. — Radial distribution of galaxies relative to the quasars in all fields which reach at least 
If =19.5. All galaxies detected at >3cr down to the average 5a magnitude limit in each field are 
included. 
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Fig. 6. — Radial distribution of galaxies in the z>lA RLQ fields. Dashed line in both figures 
is the surface density from the average published literature counts, and dotted lines are the ±lcr 
RMS uncertainties, a. All galaxies detected at >3a down to Kjjkirt = ^0.5 in the 7 z>lA RLQ 
fields with 5a limits that deep. b. All galaxies detected at >3a down to Kjjkirt = ^0 in the 15 
z>lA RLQ fields with 5a limits that deep. c. All galaxies detected at >3a down to Kjjkirt=^-^ 
in the 19 z>lA RLQ fields with 5a limits that deep. d. All galaxies detected at >3<r down to 
Kukirt=W in all 20 z>lA RLQ fields. 
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Fig. 7. — Radial distribution of galaxies in the z<lA RLQ fields. Dashed line in both figures is 
the surface density from the average published literature counts, and dotted lines are the ±lcr RMS 
uncertainties, a. All galaxies detected at >3a down to Kjjkirt=19-5 in the 5 fields with 5a limits 
that deep. b. All galaxies detected at >3a down to Kjjkirt=^ in the 8 fields with 5a limits that 
deep. c. All galaxies detected at >3a down to Kjjkirt=^&-§ in the 10 z<lA RLQ Fields fields 
with 5a limits that deep. d. All galaxies detected at >3a down to Kjjkirt=^& in all 11 fields. 
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Fig. 8. — Morphologically identified stars' isophotal r— K s colors vs. total K s magnitudes. Error 
bars are la photon noise only. Faint objects with only rightward error bars on their K s magnitudes 
have total magnitude detections of <3er, but may have isophotal magnitude detections of >3a. All 
RLQ fields are plotted except Q 1508—055 (no r data) and Q 2230+114 (nonphotometric Rc data). 
The 31 remaining fields cover 219.2 arcmin 2 and contain 380 stars detected at >3u in r or K s . 
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Fig. 9.— Galaxy r—K s colors vs. K s magnitude. All RLQ fields are plotted except 
Q 1508-055 (no r data), Q 2230+114 (nonphotometric R c data), and Q 0736-063 (uncertain 
stellar contamination). The 30 remaining fields cover 213.7 arcmin 2 and contain 3886 galaxies 
detected at >3cr in r or K s and with K s brighter than the average 5<r limits, all of which are 
plotted here. Note the reddest object at if s ~18, r—K s ^8A (3c lower limit). 
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Fig. 10. — Comparison of -fC s -selected galaxy r— K s colors between fields with actual r (or R) and 
K data (solid lines) and literature control fields with r and K data synthesized from magnitudes 
in other filters (dotted lines). Smaller histograms are galaxies with lower or upper limits to their 
colors. Each histogram has been separated normalized to unity sum. The K s magnitude range is 
given in each panel. 
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Fig. 11. — Colors of K-selected galaxies. Solid line is z=l-2 quasar fields; dotted line is control 
fields. Smaller histograms represent those galaxies with lower or upper limits to their colors. Each 
histogram has been normalized by the area imaged in the appropriate magnitude range to yield the 
surface density of galaxies per arcmin 2 in each color bin. The Kjjkirt magnitude range is given 
in each panel. The vertical scales differ between the top and bottom panels. 
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Fig. 12. — Colors of K-selected galaxies. Solid line is z=\-\A quasar fields; dotted line is control 
fields. Smaller histograms represent those galaxies with lower or upper limits to their colors. Each 
histogram has been normalized by the area imaged in the appropriate magnitude range to yield the 
surface density of galaxies per arcmin 2 in each color bin. The Kjjkirt magnitude range is given 
in each panel. The vertical scales differ between the top and bottom panels. Only a handful of 
galaxies in one quasar field contribute to the Kukirt = 2Q~21 magnitude range, so the difference 
in those two histograms is not significant. 




Fig. 13. — Colors of K-selected galaxies. Solid line is z=\A-1 quasar fields; dotted line is control 
fields. Smaller histograms represent those galaxies with lower or upper limits to their colors. Each 
histogram has been normalized by the area imaged in the appropriate magnitude range to yield the 
surface density of galaxies per arcmin 2 in each color bin. The Kjjkirt magnitude range is given 
in each panel. The vertical scales differ between the top and bottom panels. 
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Fig. 14. — Colors of If -selected galaxies in z=l-2 RLQ fields. Solid line is galaxies at 8<40" 
from the quasars; dotted line is galaxies at 9>40" from the quasars. Smaller histograms represent 
galaxies with lower or upper limits to their colors. The Kjjkirt magnitude range is given in each 
panel. 
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Fig. 15. — K/J—K color-magnitude diagram for galaxies in the control field datasets of McLeod 
etal. (1995), Elston, Eisenhardt, & Stanford (1998), and Dickinson etal. (1998). Open symbols are 
stars and filled symbols are galaxies, but classifications are only available for the McLeod dataset. 
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Fig. 16. — K/J—K color-magnitude diagram for galaxies in the five z=1.4-2 RLQ fields with J 
data. Open symbols are stars and filled symbols are galaxies. 
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Fig. 17. — J—K colors of if-selected galaxies. Solid line is 5 z=lA-2 RLQ fields with J data; 
dotted line is control fields. Smaller histograms represent those galaxies with lower or upper limits to 
their colors. Each histogram has been separated normalized to unity sum. The Kjjkirt magnitude 
range is given in each panel. 
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Fig. 18. — The Hill & Lilly (1991) richness statistic iVo.5 measured in our RLQ fields. Error 
bars are Poisson only. The horizontal dashed line at iVo.5=0 shows where the observed counts at 
<0.5 Mpc equal the prediction from data at >0.5 Mpc. 
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Fig. 19. — Histogram of "near-field" iVo.5 values measured for our RLQ fields. Abell richnesses of 
0, 1, and 2 are approximately 7V .5=10, 19, and 31. 
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Fig. 20. — The Hill & Lilly (1991) richness statistic iVo.5 measured for our RLQ fields compared to 
the literature. Error bars include the uncertainties in the corrections for fields which do not reach 
Kbcg+3. The horizontal dashed line at 7Vo.5=0 shows where the observed counts at <0.5 Mpc 
equal the prediction from the average published literature data. 
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Fig. 21. — The solid line is the excess galaxy surface density (log N) vs. Mk- Dashed lines 
show the ±1(7 uncertainty envelope. The excesses from all z<lA fields have been coadded and 
normalized by the area at each magnitude bin to make this plot. 
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Fig. 22. — The solid line is the excess galaxy surface density (log N) vs. Mk- Dashed lines 
show the ±ler uncertainty envelope. The excesses from all z>lA fields have been coadded and 
normalized by the area at each magnitude bin to make this plot. Even though these fields are at 
higher redshift, better data enables us to reach farther down the luminosity function in them than 
in the z<1.4 fields, on average. 

Fig. 23. — Color image of the field of Q 0835+580 (2=1.534) using rJK s images to drive the blue, 
green and red color guns, respectively. Saturation occurs at r=24, K s =20, and J=21 mag/arcsec 2 . 
North is up and East is left; the area in black is approximately 3.25' by 3.25'. 

Fig. 24. — Color image of the field of Q 1126+101 (2=1.516). See key to Figure ^ for details. 
The area in black is approximately 3.3' by 3.3'. 
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Fig. 25.— The SEDs for the two reddest objects in r-K s within (9=20 / /5 of Q 0835+580 (2=1.534) 
are shown as the filled and open squares with Poisson error bars on the fluxes. Data points are 
plotted at the effective wavelengths of each filter, and horizontal error bars indicate the widths of 
the filters. Model spectra are as shown on the figure and discussed in the text. 
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Fig. 26. — See key to Figure 25 for details. Model spectra are as shown on the figure and discussed 
in the text. 
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Fig. 27. — See key to Figure |25| for details. Model spectra are as shown on the figure and discussed 
in the text. 
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Fig. 28.— The SEDs for the two reddest objects in r-K s within 0=2(X'5 of Q 0835+580 0=1.534) 
are shown as the filled and open squares with Poisson error bars on the fluxes. Horizontal error 
bars indicate the widths of the filters used to construct the SEDs. The Bica etal. (1996) spectrum 
of M32 is shown as the solid line. 
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Fig. 29. — The average SED for the seven next reddest objects in r— K s within 6=2tf.'5 from 
Q 0835+580 (z=1.534) is shown as the filled squares and Poisson error bars on the fluxes. Horizontal 
error bars indicate the widths of the filters used to construct the SEDs. Model spectra are as shown 
on the figure and discussed in the text. 
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Fig. 30. — The data points are the two objects with r—K s >5 which are reddest in J—K s in the 
Q 1126+101 field, #424 (triangles) and #315 (squares), normalized to the same flux in K s . The 
lines are the 2.2 Gyr old E model of Poggianti (1997) reddened by different extinction laws with 
E(B— V)=0.35 and normalized to the data in K s . Solid line is the Milky Way (MW) extinction law 
from Seaton (1979) and Howarth (1983). Dashed line is the extinction law from Calzetti (1997). 
Dotted line is the LMC extinction law from Howarth (1983). Dash-dot line is the SMC extinction 
law given by a linear (in 1/A) fit to the data of Prevot etal. (1984). All extinction laws have been 
normalized to Ry=3. The relative extinctions for the various laws are significantly different only 
at rest-frame A < 3800 A. 
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Fig. 31. — The data points are the two objects with r—K s <5 which are reddest in J—K s in the 
Q 1126+101 field, #237 (triangles) and #18 (squares), normalized to the same flux in K s . The 
lines are the 2.2 Gyr old E model of Poggianti (1997) reddened by the same extinction laws as in 
Figure with E(B-V)=0.35. 




Fig. 32. — The data points are the two objects with J—K s >2.5 in the Q 1126+101 field which 
have believable detections in H, #425 (triangles) and #381 (squares), normalized to the same flux 
in K s . The dotted line is a 2 Gyr old GISSEL model E reddened by E(B— V)=0.75 using the 
extinction law from Calzetti (1997) and normalized to the data in K s . The dashed line is a 0.1 Gyr 
old unreddened GISSEL model E normalized to the data in r. The solid line is the sum of the two. 
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